Th ee ae ae HAWAN 
Cc 


PHILOSOPH 


— 


ICAL 
MAGAZINE 


FIRST PUBLISHED IN 1798 


«7 41 SEVENTH SERIES No. 323 December, 1950 


A Journal of 


Theoretical Experimental 


and Applied Physics 


EDITOR 
PROFESSOR N. F. MOTT, M.A., D.Sc., F.R.S. 


EDITORIAL BOARD . 

SIR LAWRENCE BRAGG, 0.B.E., M.C., M.A., D.Se., F.B.S. H 
ALLAN FERGUSON, M.A., D.Sc. 

SIR GEORGE THOMSON, M.A., v.Sc., F.R.S. 
PROFESSOR A. M. TYNDALL, C.B.E., D.Se., F.R.S. 


PRICE 10s. . 


Annual Subscription £5 2s. 6d. payable in advance, 


D AND PUBLISHED BY TAYLOR & FRANCIS LTD., RED LION COURT, FLEET ST., LONDON, E.C.4, 


Early Scientific Publication: 


DIARY OF ROBERT HOOKE, ¥.A., M.D., F.R.S. 
1672-1680 


Edited by H. W, ROBINSON and W. ADAMS 
Recommended for publication by the Royal Societ 


London 
25/- “This vivid record of the scientific, artistic and soc 
net activities of a remarkable man during remarkable years h 
too long remained in obscurity.’’—Extract from foreword 


Sir Frederick Gowland Hopkins, O.M., President of t 
Royal Society. 


MATHEMATICAL WORK OF JOHN WALLIS, D.D., F.R.S. 


By J. F. SCOTT, Ph.D., B.A. 
“His work will be indispensable to those interested int 
12/6 early history of The Royal Society. | commend to — 
students of the Seventeenth Century, whether scientific 
humane, this learned and lucid book.’’—Extract from for 
word by Prof. E. N. da C. Andrade, D.Sc., Ph.D., F.R.S. 
Recommended for publication by University of Lond 


net F 


CORRESPONDENCE AND PAPERS OF EDMOND HALLEY 


Arranged and Edited by EUGENE FAIRFIELD MACPII 


21/- First published on behalf of The History of Scien 
Society by Oxford University Press. Now re-issued 


net Taylor & Francis, Ltd. 


MEMOIRS OF SIR ISAAC NEWTON’S LIFE 


5/ By WILLIAM STUKELEY, M.D,, F.R.S., 1752 
= From an Original Manuscript 
net Now in the possession of the Royal Society, Londo 


HEVELIUS, FLAMSTEED AND HALLEY 


Three Contemporary Astronomers and their Mutual Relatic 
12/6 By EUGENE FAIRFIELD MACPIKE 
net Published by arrangement with The History of Scien 
Society 
Established 
over 150 years TAYLOR & FRANCIS, LTD 


RED LION COURT, FLEET ST., LONDON E. 
PRINTERS & PUBLISHERS OF SCIENTIFIC BOO 


CONTENTS OF No. 323. 


CVII. The Radioactivity of F2°—I. By J. V. Juttey, Cavendish 
Laboratory, Cambridge 


CVIII. The Radioactivity of F2°.—II. By R. M. Lirravmr, Cavendish 
Pe oten LO aMOTMO Ge ae eRe: fois es athe aa rs. sonienue PCRS een 


CIX. The Emission of Short Range Alpha Particles from Light 
Elements under Proton Bombardment.—IV. The Reaction 1F(pa)!60. 
By Joan M. Freeman, Cavendish Laboratory, Cambridge.............. 


CX. An Interferometric Study of some Optical Properties of 
Evaporated Silver Films. By R. C. Faust, Royal Holloway College. 
MER TEE 2. UTE... leo A ee ae aires 


CXI. The “ Velocity’’ of Brownian Movement. By R. Férru, 
iBirepecs College, University of London. 0.2. ic cenis cede cia d ewes 


CXII. The Use of Atmospherics to Study the Propagation of Very 
Long Radio Waves. By F. F. Garpner, Cavendish Laboratory, 
SWE Os ee ae eS ay aOR i ere gn ne ee 


CXIII. The Force Acting on a Body Moving Uniformly through a 
Gas containing Sound Waves. By F. R. N. Naparro, Department of 
Metallurgy, The University of Birmingham... ........ 00.2 sce seneens 


CXIV. A Note on the Decay of Current in Germanium. Diodes. By 
J. R. Titman and H. Yemm, The Post Office Engineering Research 
PELOTON Omer Pen EEN Ser oe rs Rec yartui we eos cit ateelas. Gate v's ae 


CXV. Correspondence :— 
Superconductivity of Niobium. By L. C. Jackson and. 
H. Preston-Tuomas, H. H. Wills Physical Laboratory, 
ircversiby ey DMSLOls cr tc Seen ha dae cs vs wees seems 
Lattice Sums for Ionic Crystals. By F. C. Franx, 
H. H. Wills Physical Laboratory, University of Bristol. . 
Some Extensions in Use of Resistance Thermometry in the 
Study of Gaseous Explosions. By A. 8S. Lean, 
C. RounTHWAITE and D. BrapLey, University of Leeds 


CXVI. Notices of New Books and Periodicals received :— 


J. K. Syrxix and M. E. Dyarxkina’s Structure of 
Moleculessandstne: Chemical Bond. ... 4.10. ous oes es : 


Trades they WWUITTINS oc 4a ots aon ee ernie oe or 


With Title-page, Contents, etc. 


Page 


ESS 


1214 


1281 


1292 


* * All communications for the Philosophical Magazine should be addressed, post-paid, 
to the Editors, c/o Messrs. TAYLOR AND Francis. Lrp.. Red Lion Court, 


Fleet Street, London, England. 


[ ll99 ] 


CVII. The Radioactivity of F?°.—I. 


By J. V. Jevury*, 
Cavendish Laboratory, Cambridge f. 


[Received August 23, 1950.] 


SUMMARY. 


The beta and gamma radiation emitted by F2° has been studied by 
absorption, coincidence and spectrometer techniques. The upper limit 
of the beta particle spectrum was found to be 5:03-40-05 MeV. and 
gamma rays of 1:63-0:02 MeV. and 2-45-+-0-06 MeV. with an intensity 
ratio of (8-4+2-6) to 1 were also observed. No gamma-gamma 
coincidences were found, and it is concluded that the predominant 
_ mode of disintegration of F?° leads to a 1:63 MeV. excited state of Ne29. 
The total energy release F?°—Ne?® is 6:66-0:05 MeV., and the mass of 
°° derived from this value is compared with that obtained from the 
F'® (dp) and F!® (ny) reactions. 


§ 1. INTRODUCTION. 


THE radioactive nucleus F®°, of 12 seconds half-life, decays into Ne2° 
with the emission of beta and gamma radiations. Measurements of the 
energies of these radiations and observations of beta-gamma and 
gamma-—gamma coincidences suggested the decay scheme shown in 
fig. 1. (Hornyak and Lauritsen 1948, Mattauch and Flammersfeld 1949.) 

According to this scheme the emission of a beta-spectrum with an 
upper limit of 5-1 MeV. leads to the formation of Ne? in an excited 
state of 2-2 MeV. from which gamma-ray transitions take place either 
directly to the ground state or via an intermediate state of about 1-6 MeV. 
Though the lower energy gamma-rays were known to be present (Bower 
and Burcham 1939, Curran and Strothers 1940), their energy was only 
inferred from evidence for an excited state of Ne*® at 1-6 MeV. obtained 
in observations on the neutron groups from the reaction F (dn) Ne*® 
(Bonner 1939, and Powell 1940). The total energy release in the decay 
process on the assumption of this decay scheme is in fair agreement with 
that derived from the energy release in the reaction F1° (dp) F?° (Bower 
and Burcham 1939). 

The present experiments were undertaken to improve the energy 
measurements, to obtain more precise information on the gamma-gamma 
coincidences, to determine the branching ratio for the two alternative 
de-excitation modes of the 2-2 MeV. level of Ne?® and to detect the 


* Communicated by Dr. W. E. Burcham. 
+ Now at the Atomic Energy Research Establishment, Harwell. 
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gamma-ray of energy ~0-6 MeV. which must be present according to 
this decay scheme. In the course of the work it appeared that the decay 
scheme shown in fig. 1 (decay scheme I.) was incorrect and a revised 
‘scheme is presented, together with new measurements of the beta- and 
gamma-ray energies. A new determination of the energy release in the 
reaction F1!9 (dp) F?° was also made. 

This work was completed in the summer of 1949, and a preliminary 
account of it has already appeared (Jelley 1950). It has recently been 
repeated in part and extended by R. M. Littauer in the Cavendish 
Laboratory ; his results are presented in the following paper (IL.). 


Fig. 1. 
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Accepted decay scheme of F?° (decay scheme I). 


§ 2. EXPERIMENTAL TECHNIQUE. 
(a) Production of F?°. 


F?° was produced by irradiating thick barium or calcium fluoride 
targets with deuteron beams of 1-50 wamps. at bombarding energies 
between 500 and 930 keV. The beam was interrupted periodically by 
an electromagnetically operated shutter to give suitable irradiation and 
decay periods T; and T,. During the irradiation periods, the Geiger 
counters used in measuring and monitoring the beta- and gamma- 
radiations were electrically suppressed to avoid high counting rates 
and to eliminate counting of “prompt” radiations from the target. 
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The essential features of the layout of the shutter box, analysing magnet 
and target tube have been described elsewhere (J elley and Paul 1948). 

Background counting rates, due to shutter radiations, activities present 
in the analysing magnet, and target activities of longer period than those 
‘due to F°°, were obtained by stopping the commutator used to operate 
the cyclic irradiation procedure, leaving the beam on the shutter and 
counting after a lapse of time for the F?° to decay. 

(b) Beta-ray absorption (see fig. 2 (a)). 

Two cylindrical aluminium counters, of diameter 2:2 cm., length 
6-0 cm. and wall thickness 0-034 gm./cm.? were mounted in a lead housing 
as shown, with a partition of thickness 11-6 gm./em.?. Counter M was 
used as a monitor while various absorbers were placed in front of counter A. 


Fig. 2. 
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Counter arrangements for study of F?° activity. 
(a) Beta-ray absorption. 
(6) Gamma-ray absorption. 
(c) Beta-gamma coincidences. 
(d) Gamma-gamma coincidences. 


‘The counters were situated 28 cm. from the target, T. In this experiment 
a beam of 50pamp. at 750 keV. was used, to irradiate a thick BaF, 
target, and the values for T; and T, were 10-7 and 8-8 sec. respectively. 
(c) Gamma-ray absorption (see fig. 2 (b)). 

The method of observing the secondary electron absorption (Bothe and 
Kolhorster 1929) was used to determine the energies of the gamma-rays. 
‘The F2° gamma-rays from an irradiated BaF, target emerged through 
the 0-64 gm./cm.? brass target support and entered the “ thick ” aluminium 
radiator of thickness 3-4 gm./cm.?; the 5 MeV. beta-rays were excluded 
by the radiator and target support. The secondary electrons produced in 
the radiator passed through two counters A and B in coincidence, 


A Oe 
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aluminium absorbers being placed between the counters which were of 
diameter 2-2 cm., active length 6-0 cm., and of wall thickness 34 mg./em.?. 
In this experiment T;=10-7 sec., T,;=5-2 sec., and the beam energy and — 
current were 900 keV. and 50 pamp. respectively. 

(d) Beta-gamma and gamma-gamma. coincidence experiments. 

The arrangement of counters in the beta-gamma coincidence experiments 
is shown in fig. 2 (c). A CaF, target was deposited on a 0-68 gm./cm.* 
brass support and bombarded by a collimated deuteron beam so that the 
active F2° covered an area 4X%7mm.?. To provide increased efficiency 
and to compensate for fluctuations in beam position, the gamma-rays. 
were detected by two counters connected in parallel. These counters. 
were similar to those used in the absorption experiments, but were 
surrounded by cylindrical lead converters of thickness (2-9 gm./cm.’) 
adequate to exclude the 5 MeV. beta-rays while providing the highest 
efficiency for detection of the gamma-rays. The beta-ray counter was. 
a standard cylindrical brass counter with a 10 mg./em.? mica window. 
In this experiment copper absorbers were used so that a large solid angle 
for the beta-rays could be obtained with large absorption thickness. The 
F° source activity was chosen to give suitable counting rates ; it could 
be varied by altering the beam current or energy. 

Fig. 2 (d) shows the arrangement of counters used in the gamma-—gamma 
coincidence experiment. As for the beta-gamma experiment a pair of 
counters connected in parallel were used for one channel of the coincidence 
arrangement. These (B) were two aluminium counters of the same design. 
as described above, surrounded by brass converters of wall thickness. 
3-9 gm./cm.” and mounted one each side the glass target tube. This wall 
thickness, together with the glass tube, was sufficient to eliminate the 
5 MeV. beta-rays. The other counter (A) was of special construction, 
the reasons for which will be discussed later. It consisted of another 
standard aluminium counter immersed in a graphite converter; the 
thickness of graphite between the target and counter wall was sufficient 
to exclude the 5 MeV. beta-rays. A BaF, target was used, mounted on 
a brass support of thickness 0-26 gm./cm.?, the bombarded area of the 
target was 4x5 mm.”. In both these experiments T; and T, were set 
to 1-0 and 2-6 sec. respectively. For the second part of the gamma-gamma 
coincidence experiment, in which a beta-counter A’ replaces the gamma-ray 
counter, the counter used in the beta-ray absorption experiment. was. 
employed. 

Errors resulting from variations in mean source strength during an 
experiment were minimized by maintaining the single counting rates as. 
nearly constant as possible during a run. Corrections were applied to. 
the observed coincidence rates to allow for the fact that the counting 
periods were comparable with the half-life of the activity. The 
coincidence unit was similar to that described elsewhere (Jelley 1948) ; 
the coincidence resolving time was checked from time to time and had. 
the average value 0:63-L0-02 usec. 
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(e) Spectrometer. 

The spectrometer used for the beta- and gamma-ray energy measure- 
ments was an iron-cored instrument with 90° focusing, built by R. M. 
Littauer and described in the following paper (II.). For the beta-ray 
measurements the spectrometer was connected directly to the target 
chamber in which the F?° was formed. The overall resolving power and 
geometrical aperture were 3-5 per cent and 0-031 steradians respectively 
and the instrument was calibrated for the following energies : (i) 148 keV., 
the energy of the internally converted gamma-ray F-line from ThB ; 
(ii) 6-8 MeV. which is the maximum energy of the Compton electron 
distribution from the 6-98 MeV. component of the gamma-radiation from 
the reaction F1° (p«’, y)O18 (Walker and McDaniel 1948). The calibration 

_ for gamma-ray energy measurements is described in IT. 


§ 3. RESULTS. 
(a) Beta-ray energy. 

Fig. 3 shows the absorption curve in aluminium obtained experimentally 
with the counter arrangement of fig. 2 (a). 

The curve is normalized, after correcting to a constant counting rate 
in the monitor, to an effective intensity of 100 per cent at zero absorber 
thickness. The gamma-ray background, followed out further than that 
included in the figure, was extrapolated back and intercepted the intensity 
axis at 12 per cent. The true beta-absorption curve after subtraction 
-of the background is also shown. To obtain the true origin of this curve 
a correction term, of 0-509 gm./cm.? must be added to the absorption figures 
to include absorptions in (i) the Bal, target layer, (ii) the brass target 
support, (iii) the counter wall and (iv) the air path between target and 
counters. 

The beta-ray end-point energy was derived by the method of analysis 
‘developed by Bleuler and Ziinti (1946), in which curves are constructed 
relating values of the maximum energy E,,,,, to absorber thicknesses d,, 
corresponding to reductions in intensity by factors 1/2” of that with no 
absorber ; the family of curves published by these workers strictly apply 
only to the case of allowed Fermi spectra from nuclei of Z-value 20. 
Correcting for the initial absorption mentioned above and using Bleuler 
and Ziinti’s curves, values of 4:98, 4:90, 4-92 and 4:86 were obtained for 
Esmax» giving a mean value of 4:91+0-08 MeV. A correction term of 
+1 per cent must be added in transforming from Z=20 to Z=9 and the 
corrected value for E,,,,. is then 4:96--0-08 MeV. 

The relatively small spread among the values obtained suggests that 
the beta-spectrum consists at least predominantly of one component. 

Fig. 4 shows the Kurie plot obtained with the spectrometer; the function 
4/|N(E)/f | is plotted against W, where f=1/,W?—1] and W=[1+(E/m,c’)], 
the total electron energy expressed in units of mp c*. N(E) is the relative 
counting rate per unit energy interval. From the intercept at 


W=(10-85-0-10)m,c? 
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the beta-ray end-point energy was found to be 
Bmax =5'03-L0-05 MeV. 
which agrees with that obtained by the absorption method within the 


limits of error. The beta spectrum also appears from the Kurie plot to 
consist predominantly of one component.* 


Fig. 3. 
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F°° beta-radiation absorption curve. 
Inset : more detailed survey of low absorption region. 


The period of the beta activity was found to be 12-6+-0-1 seconds 
which is rather higher than the value of 10-7+0-2 sec. recently reported. 
(Snowdon 1950). 


*In the more recent experiments reported in II., the presence of a weak 
beta component of higher energy was established. A rather higher upper limit 
(5:3 MeV.) for the main beta spectrum was obtained. 
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(b) Gamma-ray energy and intensity. 

Fig. 5 shows the secondary electron absorption curve obtained with the: 
counter arrangement shown in fig. 2 (b). The ordinates give relative values. 
of the ratio of coincidences to front counter rates (N,/N,) normalized to 
100 per cent at zero absorber thickness ; the background was obtained 
from observations with absorber thicknesses greater than those shown 
in the figure. 


Fig. 4. 
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F?° beta spectrum. Kurie plot. 


It appears from this absorption curve that there are two main com- 
ponents, with an intensity ratio (by inspection) of between 5 and 10 to 1. 
The energies of these components were derived by the method due to 
Bleuler and Ziinti (1946, fig. 8), assuming that particular intensity ratio 
(6-6/1), for the two components which gave the minimum spread in the 
energy values found for the two radiations. The energies so obtained 
were : 

E,,, 1-87, 1-75, 1-68, 1-62, 1-56, 1-51 and 1-47 MeV. 
E,,y 2°66, 2-54, 2-41, 2:36, 2-44 and 2-30 MeV. 


giving mean values 
K,4=1:640-05 MeV.., H,2=2°45+0-06 MeV. 
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It was found by experiments with gamma-ray sources of known energy 
that the ratio of coincidences to front counter rate with no absorber, 
denoted by (N,/N,)), depended on the incident gamma-ray energy. 
The ratio of intensities of the two components which gave the best 
results in the Bleuler—Ziinti analysis must therefore be corrected for the 
energy dependence of the ratio (N;/N,)y to give the true intensity ratio. 


Fig. 5. 
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F°° gamma-radiation. Absorption of secondary electrons. 
Inset : more detailed examination of low absorption region. 


With this correction, and estimating the errors in extrapolation, it was 
found that 

P= (B442-6)/1, 

L» 
which means relative intensities of 89 per cent and 11 per cent for the 
soft and hard components respectively. 
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The results obtained using the spectrometer are described in II. |The 
energy of the hard component could not be measured owing to its low 
intensity, but that of the softer radiation was found to be 1:63-L0-02 MeV., 
in good agreement with the absorption measurement. The period of the 
gamma-ray activity was found to be about 12 seconds ; no determination 
of the periods was made for the two components separately (cf. IT.). 


(c) Possible eaistence of a lower energy gamma-ray. 

According to the decay scheme L., fig. 1, there should be a gamma-ray 
of energy 800 keV. (using present values for the other gamma-ray energies) 
in cascade with the 1-64 MeV. radiation and of equal intensity. In an 
attempt to establish the existence of this radiation the low energy region 
of the curve shown in fig. 5 was studied in detail (see fig. 5 inset) using thin 
absorbers. Fig. 5 (inset) gives, together with the experimental curve, and 
fitted to it at the point for 0-15 gm./cm.?, the calculated curve for decay 
scheme I., fig. 1. Comparison of these two curves shows that there is 
no 800 keV. gamma-radiation present in intensity comparable with that 
of the 1-64 MeV. component. 

The possibility of total internal conversion of the 800 keV. radiation* 
was investigated by observing the beta-ray absorption curve carefully 
over the region 0-800 mgm./cm.? of aluminium. Thin windows were 
used on the target chamber and counters. The curve obtained is shown 
inset in fig. 3, together with the curve calculated assuming that each 
quantum of 1-64 MeV. radiation is accompanied by a 800 keV. internal 
conversion electron. The curves given by Bleuler and Ziinti were used 
for the beta-ray absorption and the data published by Marshall and Ward 
(1937) for the homogeneous conversion electrons; the calculated and 
experimental curves are fitted at the 400 mgm./cm.? point. Comparison 
of the two curves shows that there are no internal conversion electrons of 
800 keV. in the right intensity ; the slight concave rise at low absorption 
occurred with other beta-ray sources. 


(d) Beta-gamma coincidences. 

Curran and Strothers (1940) found coincidences between beta- and 
gamma-rays in the decay of F?° and and showed that the absorption of 
beta-gamma coincidences followed closely that of the beta-rays. Their 
conclusion that the dominant beta. and gamma-ray components are 
coupled, was confirmed using the counter arrangement shown in fig. 2(c). 
Observation of the single counter rates N,, N,, and the coincidence rate 
N;, were made; the latter was corrected for chance rate N, using the 
relation N,=27N,N,. A background correction, obtained by observations 
-with a thick absorber, was made to the observed values of N, to obtain 
the true counting rate in the f-counter due to F’° beta-rays; this 
correction became very large near the end-point of the beta-ray absorption 
curve. 
ee eee ee 

* The author is grateful to Professor S. Devons for suggesting this possibility. 
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In fig. 6 is shown a plot of the values of the ratio N,,/N, versus absorber 
thickness, the corrections discussed above having been applied. _ On this 
figure are also included the cut-off points corresponding to the two. 
beta-ray energies involved in the alternative decay scheme to be presented 
later. Owing to the poor statistics it was impossible to obtain evidence 
in favour of a particular decay scheme, though the result was consistent. 
with the previous evidence that at least the dominant beta ray component 
proceeds to an excited state of Ne”°. 


(e) Gamma-—gamma coincidences. 

Owing to the high energy of the beta-radiation from F*°, some or all 
of the gamma—gamma coincidences observed by Curran and Strothers. 
(1940) in their arrangement, in which lead-walled counters were used, 
may have been due to bremsstrahlung. With this in mind it was decided 
to use the counter arrangement of fig. 2 (d), in which one counter was 
constructed of graphite, in a study of gamma—gamma coincidences. 


Fig. 6. 
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Results of beta-gamma coincidence experiment. 


It is essential to carry out a gamma-gamma coincidence experiment on 
a quantitative basis if the resulting evidence is to be used to assist in the 
designation of a particular decay scheme. With this aim in view, a method 
originally proposed by Dunworth (1940) and more recently developed by 
Wiedenbeck (1947), was applied to a study of gamma—gamma coincidences 
from F°. Essentially the method consists firstly of carrying out a 
gamma—gamma coincidence experiment and, secondly, repeating it after 
one of the gamma-ray counters has been replaced by a beta-ray counter ; 
the single and beta-gamma and gamma-gamma coincidences are 
measured, from which various ratios can be set up; these ratios, which 
are different for different decay schemes, are then compared with calculated 
values, using measured values of counter efficiencies. 
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With the arrangement already described, sets of gamma-gamma and 
beta—gamma coincidence runs were carried out, counter A being replaced 
by the beta-counter A’ for the beta-gamma runs, taking care not to 
disturb the geometry of the counters B. Ten gamma—gamma and four 
beta—gamma runs were made. As before, the individual and coincidence 
rates were found from the number of irradiation cycles and counting 
periods, after a predetermined number of coincidences had accumulated ; 
these figures were set at 200 and 400 for gamma—gamma and beta-~gamma 
runs respectively. 

Correcting the single counter rates for backgrounds and the observed 
coincidence rates for chance coincidences, the following ratios were 
obtained after compounding the ten gamma—gamma and four beta-gamma 
runs :— 

From the yy runs N,,/N,,, =(3-030-23) x 10-4. \ 
From the Py runs N,,,/N,, =523+40. J 

The errors quoted here are standard deviations calculated from the 
spread among the ratios of the varoius rates in the two types of runs. 

From the above figures, the value obtained for the ratio R, defined by 


ee 
Ne Ng, 
is Rope 0°160-02. 


R is a function of the gamma-counter efficiences and the branching ratio 
of the altenative gamma-transitions ; inserting the various values for 
these quantities, the derivations of which are discussed in Appendix A, 
the value of R for decay scheme I (fig. 1) is found to be 


Ryaie= 0°47. 

It will be seen that the experiment gives too low a gamma-gamma 
coincidence rate for decay scheme I. and, furthermore, it is likely that some 
or all of these observed coincidences are false since : 

(a) an approximate calculation of the effects of bremsstrahlung showed 
that even with the graphite counter, enough coincidences could 
exist to account wholly for the observed rate ; 

(b) the formation of N® on the target, due to the reaction CP(dn)N 
could have provided a source of coincident radiations comparable 
in strength with that observed. 

A repetition of thisexperiment under better conditions (II.) has confirmed 
that the number of genuine gamma-gamma coincidences in the decay 
of F?° is very small. 

§ 4. Discussion. 
(a) The Decay Scheme of F?°. 

The small number of genuine gamma-gamma coincidences and the 
absence of a gamma-ray of energy 800 keV. in intensity comparable with 
that of the 1-64 MeV. radiation from F*° are incompatible with decay 
scheme I. (fig. 1). The results are, however, consistent with the decay 
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scheme II., fig. 7, in which the 1-64 MeV. gamma radiation follows the 
beta particle emission. The fitting of the 2-45 MeV. gamma-radiation 
into this scheme requires a complex beta particle spectrum, and a slightly 
curved Kurie plot should have been found. The expected deviation 
from linearity (shown in fig. 4) is, however, small and would have escaped 
detection. 


(b) The mass of F° 

Direction determinations of the mass of F2° in terms of that of F1* 
-can be obtained from the energy release in the reactions F1° (dp) F?° 
and F19 (ny) F2°. A repetition of the observations of Bower and Burcham 
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Proposed decay scheme of F?° (decay scheme II.). 


(1939) on the F* (dp) reaction was made (see Appendix 6), and a 
‘Q-value of 4:48-++0-04 MeV. was obtained. Using the masses of the 
‘proton, deuteron and alpha particle given by Rosenfeld (1948), and for 
the mass of F1® the value 19-00449+.0-00003 * (Tollestrup, Fowler and 
Lauritsen 1950) the mass of F®° is found to be 20-00626-L0-00005. 
The observations of Kinsey, Bartholomew and Walker (1950) on the 
energy of the capture radiation in the F!*(ny) reaction lead to a 


mass for F°° of 20-00632 +0-00006. The average of these two values is 
20-00629-+0-00005. 


ee a ae eerie aaNet ly 
* This recently published value has been adopted in preference to that used 
-in the first report of these results (Jelley 1950) 
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The present results enable the mass of F2° to be found in terms of 
that of Ne®, for which the value 19-99890+.0-00006 (Rosenfeld 1948) 
is assumed. Using the energy values E,=5-03-0-05 MeV. and. 
E,=1-63 +-0-02 MeV. the energy release in the decay of F2° is, according 
to decay scheme IT. (fig. 7) 6-66+-0-05 MeV. and the mass of F2 becomes: 
20-00606+-0-00008, which is rather lower than the previous estimates. 
If this value is correct, an additional long-range proton group might 
be expected in the F!* (dp) reaction, but this could not be found. It 
seems possible, however, from the work reported in II. that the: 
beta-particle energy may have been underestimated in the present 
experiments. 


APPENDIX A. 


Determination of the gamma-ray counter efficiencies and the evaluation of R. 
in the gamma-gamma coincidence experiment. 


Before evaluating the quantity R defined in §3 (e), it was necessary- 
to determine the gamma-ray counter efficiencies at a number of gamma-ray 
energies. We will here denote these efficiencies by « with two suffixes, 
the first indicating the energy of the radiation and the second to denote- 
‘which counters are meant (cf. fig. 2(d)). The calibration was carried. 
out by placing a 1-40 Co® source in the centre of the target area and. 
finding the single and coincidence counting rates with the counters in 
their normal position. Assuming the mean energy of the gamma-radiation: 
to be 1-2 MeV. and a zero angular correlation between the two components, 
the following absolute values for the overall efficiencies of the two counters. 
were obtained 


Use was then made of curves published recently by Maier-Leibnitz: 
(1946) to obtain the relative values of the efficiencies at other energies ; 
in this way the absolute efficiencies were determined for the following 
energies : 


Counter 0-8 16 2-4 MeV. 
A (graphite+ Al) Age (9-4 43-0 » °x 1074 
B (brass+ Al) aa 451 10S 


In this calculation it was assumed that counter A was made entirely of 
aluminium and counters B of brass; the corresponding curves given by 
Maier-Leibnitz being employed. These assumptions are justified since 
(a) the aluminium walls of the counters proper were of a thickness small 
compared to the range of the secondary electrons, (b) graphite and 
aluminium have comparable values of density and atomic number, and 
(c) over the range of energies considered the form of the efficiency-vs.- 
energy curve is similar for aluminium and brass.. 
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The ratio R, defined in $3 (e) is a function of these efficiencies and 
the relative intensities of the radiations. The value of R for decay 
scheme I. assuming no angular correlation between the various gamma-ray 
components, is given by: 

ih F(€1.68* €0-8, +£0.88’ 1-64) 
eae [f(<1-68+ €o-en) + (1 —féo-apllf(ér-6a +€o-ga) +(1 —) JeoaAl ie 
where f and (1—f) are the relative intensities of the (1-6 ; 0-8) and 2-4 MeV. 
transitions respectively. Inserting the absolute efficiencies obtained as 
described above and the value for f from § 3 (b) a value of 0-47 is obtained 
for R For decay scheme II R,,),..=0. 


R 


calc’ 


APPENDIX B. , 
Determination of the Q-value of the Reaction F1° (dp) F°. 

In this experiment a thin BaF, target was bombarded by deuterons 
and the proton groups from the above reaction detected in a proportional 
counter. Calibrated mica absorbers were used to vary the residual 
range so that the various groups could be studied separately. Another 
proportional counter with a fixed amount of absorber was used for 
monitoring purposes. Direct resolution of the groups was found possible 
using a discriminator operating at high bias. 

The protons of maximum energy were observed at three different 
angles with respect to the deuteron beam, 30°, 90° and 150° and their 
range was obtained by direct comparison with the range of the protons 
from the reactions C(dp)C'% and Be (dp) Be!® for which accurate 
Q-values are available. Throughout these experiments a deuteron 
energy of 930-+15 keV. was used. The following Q-values were obtained : 


Reference Angle of 

Reaction Observation Q MeV. 
Be (dp) Be? 30° 4-4] 
Be (dp) Be? 90° 4-4() 
Be (dp) Be!° 150° 4-44 
Cid) Ge" 30° 4-67 
C8 (dp) Ci4 90° 4-46 
C13 (dp) C4 150° 4-52 


giving a mean value of 4:48+0-04 MeV. 
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CVIII. The Radioactivity of F?°.—II. 


By R. M. LrrravEr*, 
Cavendish Laboratory, Cambridge f. 


[Received August 23, 1950.] 


SUMMARY. 


Experiments on the radioactivity of F?® are described which give the 
following results: 3-5 per cent of the disintegrations go direct to the 
ground state of Ne?°, while the remainder leave Ne*® in an excited state 
of energy 1:63-9:02 MeV. The best present estimate of the F*° mass is 
20-006 31-40-000 05, corresponding to a total energy release F?°— Ne?® of 
6-87--0-07 MeV. A weak gamma ray of energy 1-0-+0-1 MeV. is also 
observed, but its allocation to the decay of F?°is not established. Gamma— 
gamma coincidence experiments limit the fraction of disintegrations 
which produce cascade gamma rays to less than 2 per cent, unless one of 
the rays has a life-time longer than 10~ sec. 


§ 1) INTRODUCTION. 


In the preceding paper, hereafter referred to as I., J. V. Jelley has reported. 
some experiments concerned with the short-lived beta emitter F?°. 
The main conclusion to be drawn from his work is the that dominant. 
beta spectrum is linked with a single gamma ray. This result is borne: 
out in the present investigation, which was undertaken with the object: 
of clarifying some of the doubtful points raised in I. and of obtaining a. 
more precise estimate for the mass-value of F?°. 


§ 2. EXPERIMENTAL Mernop. 

(a) Production of F?°. 

°° was produced throughout this work by intermittent bombardment: 
of thin CaF, targets with 930 keV. deuterons, as described in I. The 
radiations were studied by spectrometer and coincidence counting: 
methods. In the spectrometer experiments, the irradiation and counting 
periods were approximately equal in order to achieve maximum statistical 
accuracy. In the coincidence work, where the useful source strength is 
in any case limited, the bombarding periods were shortened and a greater 
fraction of the cycle was available for counting. The period of the 
switching operations was always short compared with the life-time of 
the decay. 


* Communicated by Dr. W. E. Burcham. 
+ Now at Cornell University, Ithaca, N-Y. 
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(b) The Beta Spectrometer. 

The spectrometer is a compact instrument of the 90° magnetic focusing 
type. This is a special case of the general group discussed by Stephens 
(1934) : any homogeneous field extending over a wedge-shaped area gives 
simple refocusing of charged particles originating from a line source. 
It can easily be shown that the resolving power of such a spectrometer 
for a given transmission factor is independent of the angle of the wedge. 
The 90° instrument has the advantage over the more conventional 180° 
type that the source and detector lie outside the field gap, and that their 
separation is greater, which facilitates shielding. It has the disadvantage 
that the stray field at the edges of the pole-pieces modifies the trajectories 
of the particles. 

This effect was not considered in the design of the instrument used for 
the work on the beta particle spectrum reported in I. ; there was, moreover, 
some uncertainty in the calibration of this instrument, so that the values 
obtained with it must be regarded as only approximate. 

For the design of the present spectrometer, the stray field was mapped 
with a spinning-coil fluxmeter, and a simple graphical method was used 
to plot the trajectories of electrons on this map. In this way it was easy 
to determine the position of source and detector for optimum focusing 
in the presence of the stray field. 

Fig. 1 shows a cross-section in the median plane of the spectrometer 
(horizontal). For measurement of the beta spectra an internal target 
arrangement was used, as shown at (a). The active layer was deposited 
by evaporation in vacuo on a thin (0-5 mm.) copper backing (cooled at one 
side by water), which was mounted at 45° and had a projected width of 
4mm. The experimental resolving power with a source of this width 
is 1-8 per cent; the solid angle of acceptance of the spectrometer is 
0-6 per cent of 4. 

For monitoring the activity of F?° produced on the target, the energetic 
beta rays emerged through the backing and the window W, and were 
counted in a pair of coincidence counters not shown in the diagram. 
The spectrometer counters were also operated in coincidence to reduce 
the background due to stray gamma radiation; the shielding from the 
target was completed by building up lead blocks outside the vacuum box. 

For measuring gamma ray energies, the target arrangements shown in 
fig. 1 (b) were employed. To stop the energetic beta rays from entering 
the spectrometer, a sheet of copper (3:3 g/cm.”) was interposed. This, 

_and the thin aluminium wall of the spectrometer, acted as converter for 
Compton electrons ; no localization of this converter was possible apart 
from the fact that the bombarding beam was relatively well focused and 
thus restricted the area of material which was effective. These conditions 
were of course reproduced during the calibration of the spectrometer for 
Compton distributions. 

For an examination of photo-electrons, a thin (50 mg./cm.’) lead 
radiator, 4mm. wide, was mounted on the inside of the alumium cap. 
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The target activity was monitored by means of its gamma rays, which 
were recorded in an arrangement of coincidence counters. All readings 
were normalized to this monitor. 

The spectrometer was calibrated by reference to the following well- 
known energies: the F line of Thorium B at 148 keV., the gamma rays 
of Co® at 1/172 and 1-322 MeV. (Lind, Brown and DuMond 1949), and 
the ThC” line at 2:62 MeV. Since the iron of the magnet is far from 
saturation and the main part of the reluctance of the magnetic circuit 
lies in the air gap, the field is very accurately proportional to the magnet 
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current up to electron energies of 7 MeV. The fact that the calibration 
could be extrapolated to this region was verified by examining the gamma 
rays from F1°(px’, y)O! (at 6-13 and 6-98 MeV.; Walker and McDaniel 
1948) with the spectrometer. 

The magnet current was stabilized by an electronic stabilizer using four. 
valves of Mazda type 12E1 as control element. Relative values were 
read from a secondary standard meter, accurate to 0-2 per cent. 

(c) Coincidence Experiments. 

Fig. 2 shows the target arrangements used for the coincidence experiments 
The target was deposited on a thin molybdenum backing supported by 
two wires in the aluminium can A. (This can was lined with very thin 
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brass foil to avoid interaction of the bombarding deuterons with the 
aluminium.) The beam restricting aperture 8 was also of molybdenum, 
uncooled. Both the target and S reached dull red heat during a bom- 
barding cycle ; this prevented the formation of carbon from decomposing 
pump vapours, and thus avoided the production of N¥ by the reaction 
C™(d,n). The N' provided a troublesome background of genuine 
gamma-gamma coincidences in the earlier parts of the work, as has now 
been verified by following the decay of these coincidences (half-life 
10 mins.). 
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Target arrangements for coincidence experiments. 


Scintillation counters were used to detect the gamma rays, as they are 
more efficient than the Geiger counters used in I., and permit higher 
counting rates to be used. These counters are shown in the 180° position 
in fig. 2; to allow for a possible strong angular correlation between the 
two quanta, duplicate experiments with the counters at right angles were 
also made, but yielded no different results. An additional improvement 


4AP2 
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made in this work was the continuous recording of the chance coincidence 
rate by a method described elsewhere (Littauer 1950). This made it 
possible to use sources of such a strength that the chance rate exceeded 
the genuine rate, and also facilitated an immediate interpretation of the 
results. The coincidence resolving time was 10~° sec. 


§ 3. RESULTS. 


(a) The Beta Particle Spectrum. 

Fig. 3 shows a Kurie plot, constructed in the usual way, of the observed 
spectrum from F%°, It is obvious that the spectrum is complex, with a 
main component of end-point 5-33--0-05 MeV., and a weaker one 
(3-5 per cent) of end-point 6-740-1 MeV.* There is no evidence of 
the presence of any internal conversion line in the region 0-15—-7 MeV. 
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(A) Kurie plot from F?° beta ray spectrum. 
(B) Detailed examination of high-energy region. 


(b) The Gamma Ray Spectrum. 

The results of an examination of the Compton recoil electrons are 
exhibited in fig. 4. The experimentally determined shape due to a 
monochromatic source is shown fitted to this plot, and it is clear that the 
points in the low energy region depart significantly from this shape. 

* The energy of the main component differs appreciably from that obtained 


with the original, uncorrected gs 
i ; ed spectrometer (cf. I.) ; the present work is 
as more reliable. The earlier work did not reveal the oe energy mebcate 
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The observed differences are plotted on the same graph, and can be fitted 
reasonably well by another simple shape. It should be pointed out that 
the yields in this experiment were very low indeed and that the background, 
subtracted in Fig. 4, represented about 20 per cent of the maximum 
counting rate. The two gamma lines indicated have energies 1:64-0-03 
and 1-0-+-0-1 MeV., and intensities roughly as 4:1. Other gamma lines in 
the region 0-7-3 MeV. would have been detected if their intensity had 
exceeded about 10 per cent of that of the main component. 

An attempt to detect the photo-electrons from the main gamma 
component met with only fair success. The line appeared to be indicated 
but the intensity in this experiment was not sufficiently high to yield 
quite definite results. The gamma energy deduced from the position of 
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the line is 1-62--0-02 MeV., in good agreement with the value derived 
from the Compton distribution. The energy 1-:63-L0-02 MeV. has been 
adopted as the best estimate. 
The intensity obtainable with the spectrometer did not allow the high- 
energy gamma ray reported in I to be observed. — Better efficiency for 
gamma ray detection can be achieved with an ionization chamber using 
liquid argon, in which the Compton secondaries give up all their energy. 
G. W. Hutchinson (1948) has developed such a “liquid argon spectro- 
meter” in the Cavendish Laboratory ; at the request of the author, a 
study of the F?° gamma ray spectrum was made with this instrument. 
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The intensity was sufficiently high to allow an examination of the life- 
times of the individual gamma components, whereas in the work reported. 
in I. it was possible to check only that the average half-life of the gamma 
activity was consistent with its allocation to F”°. 

This investigation * revealed the presence of two gamma ray lines of 
energy 1-64 and 1-0 MeV. and relative intensity 4:1 approximately, in 
agreement with the spectrometer work. A complex component of energy 
up to 4 MeV., not easily resolvable into separate lines, was also found. 
The decay periods of the 1 MeV. line and of the high-energy component. 
were compared with that of the main line at 1-64 MeV., which is known to: 
belong to the F° decay (cf. I). The evidence concerning the 1 MeV. line 
was not conclusive, but did not contradict its allocation to the decay of 
F20, The period of the higher energy component, however, was definitely 
longer (~1 min.), and it is therefore likely that this component, which may 
almost certainly be identified with the 2-45 MeV. radiation reported in 
I., does not in fact arise from the decay of F?°, but may be attributed to: 
an unidentified contamination. 


(c) Gamma-gamma coincidences. 

In the analysis of the results of these experiments, we have assumed the 
efficiency of detection for gamma rays to be approximately independent 
of energy in the region 0-6—2-0 MeV., which introduces no error greater 
than about 30 per cent when the scintillation counters are used at low 
settings of discriminator bias. 

If we put 

a, 6=efficiencies in the two gamma channels, 
7=coincidence resolving time, 
nm=mean source strength, 
p=traction of disintegrations yielding two gamma rays in cascade, 
T=time of counting, 
we then have: the single counts recorded in channel A (after subtraction 
of the background) 


S(A)=2pnaT+ (1—p)nxT=(14+ p)naT. 
Similarly S(B)=(1-+-p)npT. 
The number of randon coincidences is 
R=2(1+p)2aBr . n® . T, 
and the genuine coincidences are given by 


G=2pnaBT. 
Hence 
p/(1+p)=(7/2apT)* . G/./R 
=(1/8) . G//R, 
where S=(2aBT/7)3. 


nn aera 


* The author is indebted to Mr. Hutchinson for permission to use these results 
before publication. 
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(We have assumed n?—(n)?, which is approximately true for F2° under 
. the conditions of the experiment.) 

To measure S a source of known p (e.g. Co® with p=1) was used, 
placed in the position of the target, with the same absorbers as in the 
corresponding F?° experiment. 

Table I. gives the results of the experiments as a function of the 
absorbing material used to stop the hard beta rays. 


TABLE I. 


Apparent values of » with various absorbing materials. (p is the fraction 
of disintegrations which yield cascade gamma rays). 


Absorber G/5 min. ~—_R/5 min. S p(%) 
Graphite 5+10 230 42 0-8-+1-6 
Aluminium 50-20 450 70 3-44. 1-4 
Copper 170-£40 660 120 5-8-41:3 
Lead 145+40 850 31 19+5 


It is evident that the apparent value of p rises with the atomic number 
of the absorber, and that the true value must be well below 2 per cent. 
This demonstrates the effect of bremsstrahlung from the hard beta rays, 
and also shows that there cannot be a cascade of gamma rays (of life-time 
less than 10~® sec.) in the decay of F?° with fractional intensity above 
2 per cent. 

Rough experiments were also made on beta-gamma coincidences 
from F?°. These established that there is no beta-gamma angular 
correlation greater than 7 per cent, and that there are 1-06+0-1 gamma 
quanta per beta particle. It has not been possible to separate the high- 
energy beta component for this measurement on accouut of the low 
intensity near the end-point of the spectrum. 


§ 4. Discussion. 
(a) The Mass of F?°. 

An estimate of the mass difference F29—Ne?° can be obtained (as in I.) 
by adding the energies of the prominent gamma line and of the main 
beta spectrum, for there can be no doubt that these must occur in cascade. 
Thus the total energy released in the decay is 6-96--0-06 MeV., which 
corresponds to a mass difference of 7-47+0-:06 mMU. On the basis of a 
Ne2 mass of 19-998 90--0-000 06 (Rosenfeld 1948), this gives a F?° 
mass of 20-006 370-000 08. 

Another estimate of the F2° mass can be obtained from the mass of 
Fl and the Q of either of the reactions F1%(d, p) or F1%(n, y). These 
reactions give a mean mass difference F*°—F?® of 1-001 80+0-000 04 
(cf. I.); if we take the F!® mass to be 19-004 49+-0-000 03 (Tollestrup, 
Fowler and Lauritsen 1950), then the F?° mass becomes 20:006 29+ 
0-000 05. 
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The two estimates of the F2° mass differ by 0-08--0-10 mMU, which 
indicates a satisfactory agreement. The most likely present estimate of 
the mass of F?° is 

20-006 310-000 05. 


(b) The Decay Scheme of F?°. 

If the high-energy beta component corresponds to a decay of F°° direct 
to the ground state of Ne, its maximum energy, on the basis of the 
present spectrometer measurements, should be 6-96=-0-06 MeV. The 
experimental upper limit of the distribution is 6-74+-0-1 MeV., which 
differs from the expected value by 0-22+0-12 MeV. Despite this small 
discrepancy, however, it is assumed that the transition does indeed go 
to the ground state; there is no evidence to support an alternative 
hypothesis. 

Table II. gives the ft. values (cf. Konopinski 1943) of the two partial 
beta components ; these values are not inconsistent with the assumption 
that the transitions to the excited state and ground state of Ne?® are 
respectively allowed (unfavoured) and second forbidden, as required by 
the nuclear shell model (Feenberg and Hammack 1949). On this 
assumption, and taking the excited state of Ne?® concerned to be a 
(2, +) state the ground state of F?° should be of even parity and spin 
2 or 3. 

TABLE II. 


ft. values for the two partial beta ray spectra. 

Einax t f(E) jt 
11-4 me? 12 secs. 6-5 x 10% 0-8 x 10° 
14-2 350 2x 104 O-7 >< 107 


Fig. 5 shows a number of decay schemes that might be postulated on 
the basis of these results. Scheme I, which contains only definitely 
established transitions, is incomplete in so far as it does not account for the 
presence of the 1 MeV. gamma radiation ; it must be remembered, however, 
that the allocation of this to F?° is still open to question. (The 2:45 MeV. 
component reported in I. has been omitted on the basis of the present 
results.) . 

Decay shemes IT. and III. are more speculative and would require further 
evidence to establish them. Scheme IT. postulates a beta component of 
energy about 6-0 MeV., linked with the 1 MeV. gamma ray. We have 
attempted to analyse the experimental Kurie plot into three components, 
as shown in fig. 6. The maximum intensity of the 6 MeV. component 
that can be accommodated is about 5 per cent, which is insufficient for the 
requirements of scheme II. The broken curve in fig. 6 shows the Kurie 
plot* that would be expected if scheme IT. were indeed correct. It is seen 
that little evidence to support this scheme can be adduced. 


ar 
* In the construction and analysis of the Kurie plots the effect of the finite 


spectrometer resolution has been neglected, as this is small compared with the 
other experimental uncertainties. 
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Scheme II]. (A or B) is consistent with the data from the coincidence 
experiments only if the lower of the cascade gamma transitions has a 
life-time longer than 10-6 sec., so that the quanta should not be recorded 
as coincident. If the level arrangement B is correct, so that the gamma 
transition has ouly 600 keV. of energy available, a spin >3 for the 
postulated level would suffice for this requirement. This possibility 
might repay further study. (A total internal conversion of the second 
‘quantum is ruled out by the results of both I. and the present work.) 


Fig. 5. 
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F?0 decay schemes. 
Scheme I., though possibly incomplete, is the only firmly established one. The 
energies given represent the best present estimates. 
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SUMMARY. 


The short range alpha particles emitted in the proton bombardment of 
fluorine have been examined with a magnetic analyser for proton energies 
from 300 to 950 keV. In the F(pz, y)!6O reaction three alpha particle 
groups have been measured, the average values of the energy releases 
being Q,=1-979-L0-009, Q,=1-200+0-010 and Q,—0-996-+0-010 MeV. ; 
the resonances for all three groups appear to be the same. The relative 
intensities of the groups at six resonances have been calculated. 
Angular distribution measurements indicate that there are no large 
departures from isotropy in the cases examined. Energy measurements 
of the alpha particles associated with pair emission in the F (pz, 7)!*O 
reaction give a value for the energy release of 2:056-L0-012 MeV. 


§1. INTRODUCTION. 


It is well known that the proton bombardment of fluorine can lead to the 
emission of short range alpha. particles (subsequently designated 
%,-particles) followed by gamma radiation of about 6 MeV. energy, 
according to the scheme : 


19F 41H >°Ne* >160*+4He+Q,, | 
16) * 3160) 4-y,. ; 


(1) 


Accurate measurements of the energy of the gamma rays by Walker 
and McDaniel (1948) have shown that in addition to radiation of 6-13 MeV. 
energy a component of 6-98 MeV. radiation is also present, being directly 
associated with the lower energy radiation, although the ratio of the 
intensities of the two components varies from one resonance to another. 
It seemed probable that the 7 MeV. radiation was produced in a 
competing reaction similar to (1) in which a higher excited state of 
oxygen-16 was involved : 


19f 1 17] ->20Ne* 169 ** 1 4He+ Q,, ] 
16()** 5169) 1 ,, j 


(2) 


+ Communicated by Dr. W. E. Burcham. 
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and this was proved. by observations of alpha particles of shorter range 
than those designated «, above (Burcham and Freeman 1949 a).. The 
excitation function for these shorter range particles showed resonances 
in the region of proton bombarding energies from 300 to 950 keV. 
coinciding with resonances for the production of the %,-particles, the 
intensity ratios varying considerably from one resonance to another. 
It was also found that the energy release Q, for the shorter range alpha 
particles was greater at the 873 keV. resonance than at the other 
resonances investigated. This evidence suggested the existence in the 
oxygen nucleus of excited levels at 7-15 and 6-94 MeV. above the ground 
state as well as the 6-13 MeV. level, and the pair-emitting level at about 
6-00 MeV. 

By improvement in the resolution of the magnetic analyser used for 
these investigations it has now been possible to measure the energies 
and intensities of the alpha particle groups more accurately and to show 
that at most of the resonances the shorter range alpha particles consist 
of two groups (designated «, and a). At the 873 keV. resonance the 
%-group, emitted in accordance with reaction (2), is the more intense, 
and at the other resonances studied the lower energy %-group, coming 
from 


19F} | 1H ->20Ne@* 169 ***_1 4Het Q,, | 


16 | S 
O*** 5169 14, ] 


predominates. These groups have also been studied by Chao, Tollestrup, 
Fowler and Lauritsen (1950), whose results are in substantial agreement 
with those presented in the present paper in the cases where they 
overlap}. The energy of the gamma radiation (y, and y3) associated 
with the «,- and a,-particle emission was measured carefully by 
Rasmussen, Hornyak, Lauritsen and Lauritsen (1950), and although 
the two lines were not resolved the predominant energy was shown to be 
less (6-99-+--05 MeV.) at the 873 keV. resonance than at 939 keV. where 
an energy of 7:09-+-0-06 MeV. was found. 

In the work described in subsequent sections the energies of the «,-, 
a - and «,-particles have been measured at a number of resonances 
and the corresponding Q-values deduced ; excitation functions have been 
obtained and the relative intensities of the groups at most of the 
resonances ‘from 340 to 939 keV. have been determined. Some 
observations on the angular distributions of the alpha particles are also 
presented. 

It is known that the proton bombardment of fluorine can algo lead to 
resonance reactions in which either a long range alpha particle is emitted 
or a transition may occur to the lowest excited state of 160 with the 


—— eee 


~ f We are grateful to Professor Lauritsen for communicating these-results-to 
us prior to publication. 
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emission of a short range alpha particle (designated «,) followed by an 
electron—position pair (see for example Rassmussen et al. 1950, Chao 
et al. 1950, and references included therein) 
19F 11H > 20Ne*™ + 169# 4 4He+Q 
16()*(7)_.16Q) 1 gt 4 ¢-, ; (4) 


Some measurements on the «,-particles are described in section 3(e). 


§ 2. EXPERIMENTAL MeErHop. 


The 90° magnetic analyser used for the study of the short range alpha - 
particle groups has already been described (Burcham and Freeman 
1949b, 1950a, 1950b, subsequently referred to as I., IT. and III.). In 
much of the work the magnet was used on its side so that observations 
could be made in a horizontal plane ; the proton beam which reached 
the target could be restricted to a diameter of 5 mm. and the zinc sulphide 
scintillation screen used with an EMI type 4588 photo-multiplier for 
detecting the analysed particles was then made 5 mm. wide in the 
focusing plane. Under these conditions the effective resolution of the 
system (4EH/E) was 4 per cent in energy. The analyser was calibrated. 
by observations on the peaks of proton and deuteron groups of known 
energy scattered from a thin gold film, and of deuteron and alpha particle 
groups from the proton bombardment of beryllium (cf. I.; the Q values 
used for the *Be(pd)®Be and °Be(pa)*Li reactions were 0-558 and 
2-121 MeV. respectively (Tollestrup, Fowler and Lauritsen 1949)). 

The voltage scale of the 1 MV. accelerating equipment set used for the 
experiments was calibrated by reference to the gamma ray excitation 
function for the fluorine reaction (1), observations being made with a thick 
sodium fluoride target, using a Geiger counter. The relative positions 
of the gamma ray resonances have been determined by the work of 
Bennett ef al. (1946) and Bonner and Evans (1948), and the absolute 
values were estimated by correcting their value of 862 keV. for the most 
prominent resonance in the region explored to the figure of 873-5 keV. 
found by Herb, Snowdon and Sala (1949), and multiplying their figures. 
for the other resonances by the same factor. 

The fluorine targets used for the alpha particle work were prepared by 
evaporating thin layers of calcium fluoride on copper buttons in vacuo ; 
the surface densities of the layers were estimated by weighing pieces of 
mica placed close to the buttons during the evaporation process. Target 
thickness effects were calculated on the assumption that the stopping 
power of calcium fluoride for protons of 500 keV. energy was 
260 keV./mgm./cm.? and that for alpha particles of 1 and 2 MeV. energy 
it was 1600 and 1100 keV./mgm./cm.? respectively. In most cases targets 
of thickness 0:02 mgm./cm.? or less were used. 

In the experiments on the alpha particles associated with gamma ray 
emission the yield was monitored by means of a Geiger counter, mounted 
vertically below the target ; simultaneous measurements of alpha particle 
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and gamma-ray intensities were made, the background of the Geiger 
counter being first determined by observations with a blank target. 
Counts were usually made for a standard number of microcoulombs 


collected by the target, as indicated by a current integrator. 


Big). 


Alpha particles per microcoulomb 
Alpha particles per microcoulomb 


300 400 500 600 700 800 900 1000 


Proton energy keV, 


Excitation functions for the o,-particles and the unresolved (w,-+-as)- particles 
from the °F (pa, y)'®O reaction. Left-hand ordinates refer to the region 
below 800 keV. proton bombarding energy, and right-hand ordinates to 
that above 800 keV. 


§ 3. EXPERIMENTAL RESULTS. 


(a) Alpha particles associated with gamma radiation. 

In the first experiments (Burcham and Freeman 1949 a), in which 
an analyser resolution of 10 per cent in energy was used, the lower energy 
(%- and a,-) particles appeared as a single group. The excitation functions 
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for these and for the «,-particles were observed at the same time by 
measuring the intensity at the peak of each group as a function of proton 
bombarding energy from 330 to 950 keV. Observations were made at an 
angle of 83° with the proton beam. The curves are reproduced in iated 
in which the proton energy scale has been adjusted to place the most 
prominent resonance in the «,-particle curve at the value 873-5 keV. 
The widths of the resonances are about 30 keV., which is equal to the 
estimated spread in energy of the bombarding proton beam, except in 
the case of the 598 keV. resonance in the «,-particle curve which is slightly 
broader, as expected from the observations of Burcham and Devons 
(1939), and Bonner and Evans (1948). All the resonances observed 
for the lower energy alpha particles appear to coincide exactly with the 
positions of «,-particle resonances and it may be concluded that the same 
excited levels of the compound 2°Ne* nucleus are involved in each case, 
although the branching ratios vary with proton energy. There is no 
resonance in the curve for «,- and «,-particles corresponding to the 
598 keV. resonance for «,-particles. 

For more accurate energy measurements an analyser resolution of 
4 per cent in energy was used, and the alpha particles were observed at 
an angle of 135° with the proton beam. Fig. 2(a—f) shows the alpha 
particle spectra obtained with targets of thickness ~0-02 mgm./cm.? 
at the 939, 873, 669, 598, 486 and 340 keV. resonances. Excepting the 
case of the 598 keV. resonance, where only the «,-particles appear, 
there is evidence for both «,- and «,-particle groups at each resonance, 
although for the lower resonances the position and intensity of the 
weakest («,) group cannot be fixed very accurately. During the course 
of the experiments in which these spectra were obtained a film of carbon 
contamination formed on the target surface. The development of this 
contamination was traced by observations from time to time on the 
position of the peak of the «,-particle group at the 873 keV. resonance 
and the corrections to be applied to the observed alpha particle energies 
were so obtained. The spectra shown in fig. 2 have been corrected in this 
way. Targets were frequently changed in order to avoid the necessity 
for large corrections for the carbon deposit. The effective target 
thicknesses were estimated to be equivalent to about 22 keV. for 2 MeV. 
alpha particles or 32 keV. for 1 MeV. alpha particles, so that corrections 
for target thickness were small. This was verified by observations of the 
extrapolations of the high energy sides of the spectra as described in 
TI. and IIT. 

From the corrected energies of the three sets of alpha particle groups 
the corresponding Q-values for reactions (1), (2) and (3) were obtained, 
and these are given in Table I. The probable erros are calculated from 
the estimated uncertainties in the energy calibration of the analyser and 
in the target thickness and carbon contamination corrections. The mean 
values for Q,, Q, and Q; were estimated by weighting the figures for 
each resonance according to the inverse squares of the probable errors 


Fig. 2 (a-c). 
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quoted. The values obtained (Q,=1-:979-+0-009, Q.=1-200-10:010, 
Q,=0-996-0-010 MeV.) are in very good agreement with the figures 
deduced by Chao e¢ al. (1950) from measurements at the 873 and 1353 keV. 
resonances (Q,,=1:977-0-008, Q.=1-204--0-008 and Q.,—1-002--0-008 MeV.) 
and with the value Q,=1-969 MeV. given by Strait, Van Patter and 
Buechner (1950). 


TABLE I. 


Probable Probable Q Probable 
3 


Ey Qi error Qe error error 

(keV.) (keV.) (keV.) (keV.) (keV.) (keV.) (keV.) 
939 1991 20 1199 20 999 20 
873 1980 14 1195 14 996 14 
669 1984 20 1219 25 1017 20 
598 1959 Das, — — — = 
486 1962 35 1191 3D 973 3D 
340 1979 20 1207 25 971 25 

Weighted : 
mean 1979 9 1200 10 996 10 


The relative yields at 135° of «,-, «)- and a,-particles, at each of the 
resonances studied, were obtained from the spectra shown in fig. 2 by. 
calculating for each group the integral {N(E)dE/E, where N(E) is the 
number of alpha particles counted at a given setting of the field current 
of the magnetic analyser, corresponding to the energy E. A correction 
for the number of singly-charged alpha particles expected to be 


TABLE II. 


associated with each of the observed doubly-charged groups was made 
using the data given by Briggs (1927) and Knipp and Teller (1941). 
The values thus obtained for the ratios «,/«5, «,/a3 and %4/(% +a) at each 
resonance are set out in Table II. In the cases where the ratios are 
large the figures must be regarded as rough approximations only. An 
allowance for scattered protons had to be made in estimating the yield 
of the «group at the 939 keV. resonance. In the region where the 
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results in Table II. overlap the observations of Chao et al. (1950), for an 
angle of 138°, a comparison shows general but not complete agreement. 


(b) Angular distributions. 

The angular distributions of the alpha particles associated with 
gamma-ray emission have been examined at several resonances, over 
a limited range of angles. The experimental arrangement, which has 
already been described (1, § 2 (d)), permitted observations to be made 
at angles between 60° and 135°. Thin targets were used and the variation 
in the effective target thickness as the target was rotated for each 
angular position was taken into account by the use of the gamma. ray 
monitoring counter mounted vertically below the target. Observations 
were made of the ratio of the alpha particle counting rate at the peak of 
the spectrum under examination to the simultaneous gamma ray 
counting rate (corrected for background). The distributions of «,-particles 
from the 939, 873, 669, 598 and 340 keV. resonances were observed, and 
for the energy of 669 keV. the «,-particle distribution was also studied. 
When corrected to the centre of mass system the readings at all resonances 
showed deviations from isotropy slightly outside the statistical accuracy, 
which was about 10 per cent. Now the distribution of «,-particles 
emitted at the 340 keV. resonance has been examined by Van Allen and | 
Smith (1941), and shown to be spherically symmetrical to within 2 per 
cent. It was therefore assumed that the distribution obtained in the 
present experiment at this resonance actually corresponded to an 
isotropic distribution and that apparent variations in yield with angle 
were due to systematic errors which would affect the observations at all 
resonances in the same way. The 340 keV. results were therefore used 
to correct the observations made at the other resonances ; for each case 
the yield observed at a given angle 6 was divided by the yield of 
a,-particles obtained from the 340 keV. resonance at the same angle 6, 
. The distributions derived in this way are shown in fig. 3; suitable 
normalizing factors have been applied to make the average value of the 
ratio I,(E)/I,(840 keV.) approximately unity in each case. It can be 
seen that, within the experimental accuracy (about 10 per cent) and 
over the limited range of angles which it was possible to cover, no 
significant departures from isotropy are apparent in any of the cases 

examined. This result is in agreement with conclusions of Chao et al. 
. (1950) derived on a different basis: they measured, for a number of 
resonances, the yields of «-,-, «.- and «,-particles from a given target at 
a particular angle (138°) ; the total yields over 47 steradians were then 
calculated on the assumption of isotropic distributions of all the particles, 
and these figures were found to agree to within 10 per cent with total 
yields of gamma radiation measured simultaneously. 


(c) Alpha particles associated with pair emission. . aed 
The «,-particles emitted according to the reaction (4), § I, involving 
the lave excited state of 160, were studied in the neighbourhood of the 


4Q2 
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resonance at 843 keV., which is known from measurements of the 
electron—positron pairs (Bennett et al. 1946). The alpha particle spectrum 
obtained at an angle of observation of 135° and a bombarding energy 


Fig. 3. 
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shown in the figure has an energy which would correspond to «,-particles 
from the “ tail” of the 873 keV. resonance studied above. The higher 
energy group (2-08 MeV.) is attributed to the «,-particles preceding 
pair emission according to (4). A similar spectrum was obtained at 
proton energy of 820 keV., and from the two measurements a mean 
Q-value for reaction (4) was deduced as Q,=2-056-+0-012 MeV. This 
result is in agreement with the value of 2-061-0-010 MeV. found by 
Chao et al. (1950). 


Particles per microcoulomb 


Alpha particle energy - MeV. 


Alpha particles from the reactions *F(pa, y)'®O and F(pa, z)'®O at a proton 
bombarding energy of 830 keV. 


§ 4. DISCUSSION. 


From the present state of knowledge concerning reactions produced. 
by the proton bombardment of fluorine one can draw up a generalized 
disintegration scheme as shown in fig. 5. The excited states formed in 
20Ne appear to fall into two classes : (i) levels from which transitions may 
occur either to the ground state of 1°O with long range alpha particle 
(a9) emission or to the lowest excited state of °O with short range alpha 
particle («,) emission followed by the emission of an electron—positron 
pair (7) ; (ii) levels from which transitions may occur to any one of the 
next three excited levels of 1®O with emission of short range alpha particles 
(a, % or a) followed by gamma radiation (yj, y2 or ys). In the latter case 
the relative probabilities of the transitions to the three 1°O levels vary 
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‘considerably from one resonance to another as is shown by the intensity 
‘ratios listed in Table IJ. An attempt was made to determine the way 
in which these relative probabilities might be governed by the 
penetrability factors for the emitted alpha particles, whose energies are 
well below the barrier of the residual #*0 nucleus. The penetrability 
‘factors were calculated from the formule given by Bethe (1937), assuming 
a constant value for the radius of the 60 nucleus derived from the data 
of Amaldi and Cacciapuoti (1947). The factors appropriate to alpha 
particles emitted with angular momenta of all values up to l=4 were 
calculated and an attempt was made to select consistent sets of the ratios 
of these factors which would fit the observed «/«, and «/a, ratios at the 
various resonances, and give some information about the angular 


Fig. 5. 
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momenta of the emitted particles and the spins of the 160* levels. No 
such self-consistent sets could however be found. One is tempted to 
conclude, in a general way, that in the cases where the ratios « fxg or 
a,/4; are fairly small the angular momenta of the «,-particles are erent 
than those of the lower energy particles, so that the emission probabilities 
of the latter are raised to comparable values. On the other hand. in the 
case of the 598 keV. resonance the high relative probability of « particle 
emission could be accounted for if these particles were ane to be 
emitted with small or zero angular momentum while the lower energy 
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particles encountered an appreciable angular momentum barrier. Such 
an assumption would also account for the greater width of the 598 keV. 


resonance compared with the other resonances (cf. Devons and Hine 
1949). 


The recent angular correlation measurements made in this laboratory 
for the «,-particles and y,-radiation from the 340 keV. resonance (Barnes, 
French and Devons 1950) indicate a useful method of obtaining data 
on the angular momenta of the «-particles from the various 2°Ne* levels, 
and the spins of these and the 1%O0* levels. When such data have been 
obtained, further deductions may then be possible from the observed 
intensity ratios discussed above. It may appear, for example, that the 
radius of the 1°O* nucleus is not the same for the three levels of 
excitation involved. 
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ABSTRACT. 


Thermally evaporated silver films have been studied in the reflectivity 
range of 4 to 95 per cent. The air-silver and glass—silver reflection and 
transmission coefficients, the corresponding reflection phase changes and 
the apparent film thickness have been measured and compared with those 
of sputtered silver films. For very thin films, in a region where the 
light absorption is high, the glass-silver reflection coefficient passes 
through a minimum, and in the yellow and the green, but not in the blue, 
the glass-silver reflection phase change suddenly alters from an advance 
to a retardation. This transition in the phase change is intimately 
associated with the reflected fringe asymmetry. 


§1. IyTRODUCTION. 


‘Tas presence of metallic coatings on the plates of an interferometer 
modifies the fringe patterns, the reflected system being affected differently 
from the transmitted one. In fig. 1 are seen three isotropic dielectric 
media 1, 2, 3, separated by the two parallel planes XX’ and YY’ which 
may or may not be coated with metallic films. The central medium has a 
refractive index and a geometrical thickness t. Monochromatic radiation ~ 
which is incident upon XX’ enters medium 2 at an angle #@ with the normal, 
and gives rise to a transmitted series of parallel wave trains a, b,c, ..., 
and a reflected series A, B,C,.... The optical properties of the boundary 
XX’ for light incident from medium 1 are an amplitude reflection 
coefficient 7,, an amplitude transmission coefficient ¢,, a reflection phase 
change «, and a transmission phase change 8, ; if the light is incident 
upon XX’ from medium 2 the same symbols are used with the suffix 2. 
If the light falls upon YY’ from medium 2 the reflection and transmission 
coefficients are p and +, whilst the reflection phase change is y. These 
optical values, which apply to light of vacuum wavelength A incident 
upon the interferometer at one particular angle, are in general dependent 
upon the state of polarization. 

a 

* Communicated by Professor S. Tolansky. 
+ Now at the British Rayon Research Association, Urmston. 
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The resultant intensity of the transmitted light is a periodic and 
symmetrical function of the phase quantity 


d=4n/A.ptcosO+a,+y . . . . . . . (1) 
the greatest and least intensities occurring at 5=2na and (2n-+ 1) 
Tespectively. It has, however, been shown by Hamy (1906) that the 
reflected intensity is a periodic but generally asymmetrical function of 8. 
A given reflection minimum is therefore not necessarily equidistant 
(in terms of 5) from the two neighbouring reflection maxima, nor does 
the condition for a reflection minimum generally correspond to that 
leading to a transmission maximum. It is to be noted, however, that 
whenever XX’ is an uncoated dielectric interface the reflected fringe 
pattern is symmetrical, the reflected intensity minima coinciding with the 
transmitted maxima. 


Bigs 1: 


Hamy confirmed his theoretical predictions by observing the reflected 
fringe system produced by interferometers whose plates had been coated 
with chemically deposited silver films. More recently silver films 
evaporated on to glass have been similarly studied by Holden (1949), who, 
having compared his findings with the optical properties of sputtered 
films as determined by Rouard (1937), suggested that certain features of 
the reflection fringe asymmetry are closely associated with the glass-silver 
reflection phase change. In order to test this idea the properties of 
evaporated silver films have been examined at normal incidence and 
compared with those of sputtered films. The apparent film thickness, 
the reflection and transmission coefficients and the refiection phase 
changes at both the air-silver and glass-silver interfaces were measured 
and the results related to the observed reflection fringe asymmetry. 

The silver, evaporated from molybdenum strip filaments at a pressure 
of about 510-5 mm. Hg, was deposited on to two glass optical flats 
and a thin glass eoverslip. The variation in the mass deposited upon 
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unit surface area was so slight that the three silverings could be regarded 
as identical. The two flats formed an air-gap interferometer, whilst the 
coverslip was employed for the determination of the glass-silver reflection 
phase change. All the measurements were made in the space of a few 
hours during which time no significant changes in the films’ properties 
or appearance were detected. 

The electron microscope has revealed a discontinuous structure in 
semi-transparent metal films, and thus one may not strictly speak of the 
thickness of a film. The interferometric method |Khamsavi and 
Donaldson (1947), Avery (1949)] which was employed yields an apparent 
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thickness. D, which is that the film would have if it were uniformly 
distributed and of the same density as the bulk metal. A strip of the 
evaporated film was gently removed with a razor blade and an opaque 
silver film deposited over the whole surface. the height of the well- 
defined step being determined from reflected Fizeau fringes. The 
probable error in the mean step height did not exceed +15 A. for apparent 
film thicknesses in excess of 100 A. and +25 A. for thinner films. 
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§2. THE REFLECTION AND TRANSMISSION COEFFICIENTS. 


An electron multiplier cell was used to compare the incident light 
intensity with the intensities of the beams reflected and transmitted by 
the silver films, the angle of incidence being only 4°. After the 
introduction of correction terms necessitated by the presence of the 
glass substitute, the desired air-silver and glass-silver reflection and 
transmission coefficients were calculated. Three wavelength bands were 
selected by filters; 6700 A. to 5900 A. with a peak at 6300 A., 5500 A. 
to 5200 A. with a peak at 5400 A., and 4650 A. to 4350 A. with a peak at 
4500 A. 


100 200 300 400 D—> 500 


The results are presented graphically in fig. 2. As the air-silver and 
glass-silver transmission coefficients agreed to within the estimated error 
(+0-02 T), a mean value T has been taken. The error in the air-silver 
reflection coefficient R is about -+0-01 R for high reflectivities, but 
increases to +0-04 R when the reflectivity is below 10 per cent. The 
inaccuracies in the glass-silver reflection coefficient R’ are greater, rising 
from +0:02 R’ to as much as +0-5R’ at the reflectivity minimum. 
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The absorption for light incident from the air side A and from the glass 
side A’ have been calculated from the measured quantities, the presence 
of any light scattering being neglected. Being difference terms their 
errors are relatively high and small fluctuations in the calculated absorp- 
i are not necessarily significant. ; 
Tees results have aly anes with those obtained by Goos (1936) 
and Rouard (1937) for sputtered films, and by Krautkramer (1938), 
Strong and Dibble (1940), and Sennett and Scott (1950) for evaporated 
silver films. The following apparent general conclusions may be drawn :— 
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(1) At a given wavelength the air—silver reflection coefficient R decreases 
with diminishing film thickness until the value appropriate to the substrate 
alone is reached. 

(2) The wavelength variation of R depends upon the film thickness. 
Normally R increases with increasing wavelength, but with very thin 
films the opposite behaviour has been found by Goos, Sennett and Scott, 
and the writer, whose 25 A. thick film had the coefficients 5-0 per cent. 
{6500 A.), 6-9 per cent. (5400 A.) and 7-4 per cent (4500 A.). 
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(3) At a given wavelength the substrate-silver reflection coefficient R’ 
decreases with diminishing film thickness until a low minimum value is 
reached ; it then increases finally attaining the value appropriate to the 
uncoated substrate. R’ is apparently always less than R. 

(4) For thick films R’ rises with increasing wavelength, but for thin 
films in the neighbourhood of the R’ minimum the reverse may occur, 
this behaviour being responsible for the beautiful reflection colours 
observable in thin films. The 50 A. film, for example, was blue-mauve 
in reflection from the glass side, R’ being greater in the blue than in the 
red. The 25 A. film appeared rose-pink, R’ now being weaker in the blue ; 
the reflection colour from the air side of this film was a brilliant torquoise 
because R possessed the relatively uncommon property of being strongest 
in the blue. 

(5) The transmission measurements made by the different workers 
show little agreement, the major differences occurring in the neighbour- 
hood of the substrate-silver reflectivity minimum. In this region the 
writer has found that the general increase in transmission with 
diminishing film thickness is interrupted by a marked local minimum. 
This was not observed by Goos or by Strong and Dibble. (Rouard and 
Krautkramer did not take transmission measurements in this region.) 
The recent work of Sennett and Scott is very illuminating ; they found a 
minimum with high but not low rates of evaporation, the relevant factor 
being the film structure. 

(6) Whether there is a minimum in the transmission curve or not, 
the absorption terms A and A’ have peak values for all wavelengths in 
the neighbourhood of the substrate-silver reflectivity minimum. The 
peak absorption for both A and A’ decreases with shortening wavelength, 
but for the purpose of clarity A’, which is everywhere greater than A, 
has not been plotted. According to Sennett and Scott the slower the 
rate of absorption, the greater the film thickness for which the peak 
occurs, the onset of this high absorption being closely related to the 
state of aggregation of the film. For evaporation times of less than a 
minute the maximum absorption was found for an approximate film 
thickness of 100 A. ; the figure obtained by the writer using high evapora- 
tion rates was 60 A. (The slight dip near the peak of the absorption 
curve lies within the experimental errors.) 

Tt is thus evident that the exact mode of film preparation plays a very 
important role in determining the film properties, and quantitative 
comparisons between sets of films made under different conditions are 
unsatisfactory. The general qualitative behaviour described above is, 
however, apparently common to all sets of silver films, evaporated or 
sputtered. 

§3. OBSERVATION OF THE FRINGE PATTERNS. 

The theory of reflection fringe asymmetry was originally developed 
for a plane-parallel sheet interferometer, but both v. d. Pahlen (1912) 
and Brossel (1947) have shown that provided the gap thickness and the 
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slope of a wedge interferometer are small, the fringe intensity distribution 
on a Feussner surface differs negligibly from that produced by a sheet 
interferometer. A wedge can therefore be utilized to study the fringe 
patterns produced by metal-coated interferometer plates, any asymmetry 
due to the wedge slope being insignificant. There are two zero order 
Feussner surfaces on which the fringes lie, one in transmission, the other 
in reflection, but they do not exactly coincide even when the light falls 
upon the interferometer at normal incidence. In order therefore that 
both surfaces could be focused simultaneously a well-collimated normally 
incident beam, an interferometer of small gap and slope, and a 2-inch 
(N.A. 0-17) objective were used. 

Fizeau fringes and fringes of equal chromatic order (Tolansky 1945) 
were studied in both transmission and reflection. The observation of 
the fringes of equal chromatic order was carried out, using a specially 
constructed plane grating spectrograph (14,200 lines/inch, inverse dis- 
persion 50 A./mm.). The spectrum was normal, wavelength differences 
being proportional to the linear displacements on the photographic plate. 


§4. THe Atr-SILVER REFLECTION PHASE CHANGE. 


Consider an air-gap interferometer with an unsilvered front surface XX’ 
and a rear surface YY’, of which one half is silvered. The reflection 
Fizeau fringes, the minima of which coincide with the transmission 
maxima, run from the glass/glass to the glass /silverregion and at the junction 
a fringe displacement is observed (figs. 11, 13, Pl. XX XVI.). At normal 
incidence the positions of the reflection minima in the two regions are 
determined by the equation 


N= 21 OAS Oe a i 


in which the dash notation applies to the unsilvered region. This is 
derived from equation (1), t and ¢’ being the two gap thicknesses at 
which a fringe of given order n occurs. Q and Q’, representing the 
phase terms, are related by (Q—’)= 1/27. (y—y’), where y and y’ are the 
phase changes (in radians) at the air-silver and the air—glass surfaces 
respectively. 

If dqA represents twice the silver film thickness and dp the observed 
fractional order shift of the two fringe systems, then one obtains 


2(t —t)=(dq+dp)A=A2n.(y—y’), . . . . . (3) 


in which dp is positive when the fringes in the silvered region are on the 
side of decreasing gap thickness relative to those in the unsilvered region. 
It is known from the Fresnel coefficients that the reflection phase change y’ 
at an air-dielectric interface is +7, the sign being in doubt. Following 
Rouard (1937) the negative sign is here accepted, the phase change being 
interpreted as an advance, 7.e. the optical gap is smaller than the 
geometrical one. 
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If a film is of high reflectivity the reflection fringes are barely visible 
in the region where the silver film covers the rear surface. Transmission 
fringes can be used instead, the fringe contrast and sharpness being 
enhanced by coating the front surface with a highly reflecting, semi- 
transparent film. Equation (3) also holds for these transmission fringes 
which are unfortunately very diffuse (cf. fig. 3, Pl. XXXVI). 

Such measurements were made for the three mercury radiations 
5780, 5461 and 4358 A. The phase changes, all of which are advances, 
are shown graphically in fig. 4 as fractions of 7. For film thicknesses 


Fig. 4. 
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above 100 A. the error is -+0-02, whilst for thinner films it reaches 
+0-03. With a given film the relative values for different wavelengths 
are more accurate (+0-01), since they are uninfluenced by uncertainties 
in the film thickness. 

An alternative method of determining the air-silver phase change is 
to study the transmission fringes of equal chromatic order arising within 
an air-gap interferometer formed from two identically silvered optical 
flats Such a fringe pattern is shown in fig. 5 (Pl. XXXVI.) for 529A 
thick films, the mercury spectrum being superimposed upon one half of 
the field of view. The wave-number difference between any two trans- 
mitted fringes of order n and (n-+1) is 


Nieenn 1 
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the order increasing with decreasing wavelength. As the variation of 
Q from order to order (i. e. with wavelength) is slight, the wave-number 
difference closely approaches the constant value }t. The wavelength 
difference (A,,—A,,,,) must therefore decrease rapidly as the wavelength 
shortens. The effective order of interference may be found from the 
equation 
r 
n—Q=2t/A,= —“SF- ww ws se (5h 
Na Nad 

in which, as N. Barakat has remarked, the fractional part represents 2, 
n being of necessity integral. Here 2=y/7 where y is the desired air— 


silver reflection phase change. 
The fringe pattern of fig. 5 (Pl. XX XVI.) is analysed in Table I. The 


TaBeE [. 
Wave- Wave- Wave- Wave- Q 
length number length number 6.4v n—Q Fi 
(A.U.) (om -*) differences _ differences eas 
6457 15487 270 676 — 22-83 — 
6187 16163 247 672 4037 23-82 — 
5940 16835 229 675 — 24-82 -= 
5711 17510 211 672 4035 25-81 —0-81 
5500 18182 196 672 a 26-80 —0-79 
5304 18854 181 670 4033 27-79 — 
5123 19524 LUE 674 _- 28-77 a 
4952 20198 160 670 4028 29-77 - 
4792 20868 149 670 — 30-76 os 
4643 21538 139 664 4020 31-75 — 
4504 22202 132 671 + 32-72 — 
4372 22873 125 673 4019 33°71 —0-72 
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wave-number differences are almost constant but, as seen from the fifth 
column, there is a gradual diminution due to the slowly varying phase 
change. The geometrical gap thickness can only be found from the 
approximate equation 4dv—4jt, and as a small error in ¢ leads to a large 
percentage error in 2, another procedure had to be adopted. The value 
chosen for ¢ was such that the phase changes calculated from the fringes of 
equal chromatic order agreed as closely as possible with those determined 
from the Fizeau fringes. In this particular example the chosen value 
was 2t—1-474 x 10-3 cm. (dv=678 cm.~1), the calculated effective order 
being given in the sixth column. In the final column is the value of 
{2 obtained from Fizeau fringes. To ascribe a unique value to the 
geometrical gap implies an abrupt interface between the silver film and 
the air, and is therefore equivalent to defining the film thickness. In 
this way the phase changes at the two wavelengths 6300 A. and 4000 A. 
were found for films thicker than 125 A. These are represented by the 
dotted curves in fig. 4. 
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An opaque film, whose phase change was required, was matched 
against a semi-transparent 370 A. film and the reflection fringes of equal 
chromatic order were studied (fig. 6, Pl. XX XVI.). It was known that the 
reflection minima coincided with the theoretical positions of the trans- 
mission maxima (cf. §6), and the phase change could therefore be 
calculated from equations (4) and (5) in which it was possible tointroduce an 
approximate value for the correction term (2, ,,—Q,). Here, however, 
Q=1/27 . (x.+y) where «, is the already-determined phase change at the 
transparent film and y is the desired phase change at the opaque film. 
The values obtained are close to those found for relatively thick 
semi-transparent films. 


Wavelength (A.) 6300 5780 5460 4360 4000 
y (in terms of 7) 0-84 0-82 0-80 0-73 0-70 


Several conclusions may now be drawn from fig. 4. 

(1) For thick films the air-silver reflection phase change increasingly 
departs from 7 with shortening wavelength, the difference in the phase 
change between 6300 A. and 5461 A. being almost the same as that 
existing between 4358 and 4000 A. These effects, corresponding to the 
more rapid falling off of the reflectivity as the wavelength is lowered, 
were also detected by Meggers (1921) for cathodically sputtered, thick 
‘silver films. 

(2) As the film thickness is reduced the phase change oscillates, 
a well-marked drop occurring at about 125 A., 7. e. for a reflectivity in 
green of about 25 per cent. Such a minimum was also noted by Rouard 
for a film thickness of 90 A. and a reflectivity of about 20 per cent. In 
this region the wavelength sensitivity of the phase change is very slight. 

(3) For very thin films the phase change is close to 7 and even becomes 
greater than this quantity, thus being observed as an advance relative 
to that at a dielectric face. Such advances were also measured by Rouard 
for films below 32 A. thickness. With such thin films the phase change 
deviates least from 7 in the blue in accordance with the unusual behaviour 
of the air-silver reflection coefficient. Further reductions in the film 
thickness cause the phase change to attain the value 7 appropriate to 
the substrate alone. 


§5. THE GLASS-SILVER REFLECTION PHASE CHANGE. 

The glass-silver reflection phase change was measured by viewing 
reflection Fizeau fringes within a thin glass cover slip, the rear face of 
which was half-covered by the silver film. The phase change at the 
glass-air face, which according to theory is zero, was accepted as a 
standard. A knowledge of the film thickness was not required in order to 
calculate the glass-silver phase change which is provided by 27 .dp, where 
dp is the observed fringe displacement at the boundary of the two halves 
of the interferometer (cf. equation (3)). When the film reflectivity 
was too great to yield satisfactory reflection fringes the front surface of 
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the cover slip was coated with a highly reflecting, semi-transparent film 
and the transmission fringes were then studied. The lack of planeness. 
of the faces of the glass slip often produced an irregular fringe spacing, 
and the very low glass-silver reflectivity of thin films resulted in a poor 
fringe visibility. For these reasons an accuracy better than +-0-04 cannot 
be claimed for the phase change measurements shown in fig. 7 as 
fractions of 7. 

For thick films the phase change is an advance which more closely 
approaches 7 the longer the wavelength. In the yellow and green a 
decrease in film thickness causes this advance to become a marked 
retardation, the transition, which covers a very narrow thickness range, 


Hig. 7, 


gra 


occurring at that thickness (60 A.) for which the glass-silver reflection 
coefficient is at its minimum value. This retardation then diminishes, 
presumably reaching zero, the value appropriate to the uncoated 
dielectric. In the blue, on the other hand, the initial phase advance 
decreases and might become a very slight retardation before it finally 
vanishes. Similar results were found with the sputtered films studied by 
Rouard, the conversion of an advance into a retardation occurring in 
the neighbourhood of the glass-silver reflectivity minimum for ae 
lengths greater than 4358 A. 

The reflection fringes obtained with the 61 A. film at the wavelength 
5461 A. are shown in fig. 8 (Pl. XXXVI.), the fringes in the silvered part 
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(the lower region of poorer contrast) being on the side of increasing gap 
thickness relative to those in the unsilvered part. The glass-silver phase 
change is therefore an advance. The film edge is not abrupt, and as the 
deposit becomes thinner this advance rapidly diminishes and even becomes 
a retardation, the fringes moving to the side of decreasing gap. This 
demonstrates the high sensitivity of the phase change to alterations in 
the thickness of such low reflectivity films. 

The glass-silver phase change apparently deviates from 7 to a greater 
extent than the air-silver phase change for all wavelengths and film 
thicknesses, this behaviour probably corresponding to the relatively 
low reflectivity on the substrate side of the film. 


§6. REFLECTION FRINGE ASYMMETRY. 


The phase quantity y, which is determined solely by the optical 
properties of the upper interferometer surface XX’ according to the 
expression 


Vag 8 ee ey ee ee (0) 


decides whether or not the reflection fringes shall be asymmetrical. 
It is, however, the combined properties of both interferometer plates which 
shapes the reflection fringe pattern and determines at what value of 
5 the reflection minima will occur. The theoretically predicted fringe 
behaviour is summarized in the following Table, in which it is to be noted 
that =2n7 always corresponds to the transmission maxima and that 


[= 2tytep[p?ro(7yrot+tyte)—1 1]. 


TABLE II. 
cys ; ‘ Reflection minima 
x Condition Reflection fringe pattern Se ates 
(2N-+1)z <0 Symmetrical 2na 
(2N+1)z A) Uniform field oa 
(2N-+1)z I>0 Symmetrical (2n-++1) 
2Na Symmetrical (2n+1)a 
Between (2N+1)z Asymmetrical, minima closer Between 2n7 
and 2Nz to maxima of lower order and (2n-+ 1): 
Between (2N+1)z Asymmetrical, minima closer Between 2n7 
and (2N+-2)z to maxima of higher order and (2n—1)z 


For relatively thick films the form of the reflection fringes was studied 
by simultaneously viewing the reflected and transmitted fringe systems. 
In this way the Fizeau fringes of fig. 9 (Pl. XX XVI.) were obtained with the 
529 A. thick films at the wavelength 5461 A., faint secondary fringes arising 
from the 45° reflector being visible across the reflected pattern. The 
reflection minima agreed in position with the transmission maxima to within 
0-02 of an order, a difference attributable to collimation and focusing 
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uncertainties. Fringes of equal chromatic order were also observed with 
the same interferometer, the close agreement of the fringe positions 
throughout the visible spectrum being seen from fig. 10 (Pl. XXXVI.). This 
coincidence of the reflection minima with the transmission maxima was also 
discovered for all film thicknesses in excess of 142 A. (47 per cent reflectivity 
in the green), nor was any asymmetry perceptible. The phase term x 
must, therefore, be close to an odd multiple of z for all wavelengths. 

Because of the diffuseness of transmission fringes at low reflectivities 
a different method was adopted for thinner films. Only a half of each of 
the two glass interferometer plates was silvered and the interferometer 
was divided into four regions by the intersection of the two straight line 
edges of thefilms. In fig. 11 (Pl. XX XVII.) are seen the four separate Fizeau 
fringe systems observed in reflection (5461 A.) with the 93 A. films. Quoting 
the upper interferometer face first the four regions are glass/glass, 
glass/silver, silver/glass and silver/silver ; the reflection minima in the 
two first-named regions coincide with the transmission maxima at 6=2n7 
and thus form reference marks. The fringe displacements for the 
93 A. films are given below as fractions of an order, the fringes in the 
second region being on the side of decreasing gap thickness, 7. e. 
decreasing 6. 


A.  glass/glass—glass/silver 0-067 
B.  glass/glass—silver/glass 0-147 
C. _ silver/glass—silver/silver 0-067 
D.  gilass/silver—silver/silver 0-152 
E.  glass/silver—silver/glass 0-080, 0-085. 


The origin of the displacement A as the combined effect of the film 
thickness and the air-silver reflection phase change has already been 
described in §3. The inequality of the displacements A and B 
shows that the addition of a silver film to the upper plate XX’ is not 
equivalent to adding one to the lower plate YY’, the difference 
EK representing the influence of y. A value for E can also be obtained 
from displacements C and D. As, however, the positions of the reflection 
minima are dependent upon the reflectivity of the lower interferometer 
face whenever y deviates from a multiple of 7, the displacement A might 
differ from C and displacement B from D ; the two calculated values of 
E should, however, be identical. No conclusive evidence was found for 
such differences, this probably being due to the small alteration in the 
rear surface reflectivity (4 to 14 per cent) from one region to another. 

In fig. 12 are seen the differences E expressed as fractions of an order, 
positive values signifying that x is larger than (2N-+-1)z and negative 
ones that it is below (2N-+-1)z. When E is +0°5, x is an even multiple 
of 7; this condition is closely approached with the 50 A. films in the 
green, the reflection fringes in the silver/glass and silver/silver regions 
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then possessing the transmission-like appearance illustrated in fig 13 
(Pl. XXXVIT.). When y is no longer an integral value of 7 the reflected 
fringe pattern is asymmetrical in the sense indicated in Table IL. 

‘The 61 A. thick film is of interest In the yellow and the green the 
Fizeau fringe asymmetry is such that the minima are closer to the 
maxima of higher order. The displacement E, however, is apparently 
negative, not positive, since the minima in the silver/glass region appear 
to be on the side of greater order relative to those in the glass/silver 
region. ‘Two explanations are possible for this behaviour : y lies between 
(2N-+-1)z and (2N+2)z and some spurious edge effect gives the impression 
of a fringe movement towards larger orders instead of smaller ones, or 
x is between 2Nz and (2N—1)z, in which instance the displacements 


Fig. 12. 


X =(2N+t2)7---- -- 


X= (2N+1)7 


E in the yellow and green should read —0-58 and —0-69 instead of 0-42 and 
0-31 The first explanation is favoured, the spurious influence being pro- 
vided by the rapid variations in the glass-silver phase change which take 
place at the film edge (cf. fig. 8, Pl. XX XVI.), these variations bringing 
about marked changes in the value of y in accordance with the arguments 
of the following section. In the blue, where the glass-silver phase change 
suffers only mild alterations and does not change sign, no edge effects 
are noticeable, and E is undoubtedly positive. This strong wavelength 
sensitivity of the film properties is illustrated in fig. 14 (Pl. XX XVIL.), 
which shows fringes of equal chromatic order taken across the boundary 
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between the silver/silver and the glass/silver regions of an air-gap 
interferometer formed with the 61 A. films. In the blue the silver/silver 
minima are close to but on the low order (long wavelength) side of the 
glass/silver ones. As the wavelength increases the silver/silver minima 
become more displaced relative to those in the other region, but despite 
appearances this fringe displacement is towards longer and not shorter 
wavelengths, the strong slope of the silver/silver fringes at long wave- 
lengths representing marked changes in the optical properties as the film 
thickness alters slightly. 


$7. DISCUSSION. 


If the film thickness is greater than about 150 A. the phase term x is 
close to (2N-++1)z for ‘all wavelengths, this being the value appropriate 
to the dielectric substrate. At the wavelength 4358 A. further reductions 
in film thickness cause a small initial increase in y which is then followed 
by a decline, (2N-+1)z again being reached for the thinnest films. The 
behaviour at the two longer wavelengths is, however, different. At 
first y increases rapidly, closely approaching (2N+2)z in the yellow and 
possibly also in the green. Then at an approximate thickness of 60 A. it 
abruptly drops to about 2Nz ; this rapid decline is followed by an increase, 
(2N-+-1)7 being reached from the lower side and not the upper one as in 
the blue. Holden (1949) who has utilized microphotometer traces to 
study the reflection fringe asymmetry in the green, discovered that the 
magnitude of x rapidly diminishes from about (2N-+2)7 to 2Nz at an 
approximate air-silver reflectivity of 15 per cent, this being in agreement 
with the present work where the transition occurs at a reflectivity in the 
green of 12 per cent. 

In equation (6) the quantity y was defined as the sum of the air silver 
and glass-silver transmission phase changes minus the sum of the air— 
silver and glass-silver reflection phase changes. The two transmission 
phase changes are unknown and their contribution cannot be assessed. 
It is reasonable to assume that they suffer no rapid alteration in magnitude 
within the thickness range 50 to 70 A. In the same region the air-silver 
reflection phase change undergoes relatively mild oscillations, and the 
glass-silver reflection phase change must therefore exert the greatest 
influence. For example, at 60 A. this phase change abruptly increases 
by about 1-47 in both the yellow and the green (phase advances are 
negative in sign). At the same thickness the fringe displacement 
E decreases by about 0-75 of an order, thereby indicating that y has been 
reduced by an amount between 7 and 27. In the blue, however, both 
the phase change and y show only gentle variations and neither reverses 
sign. The general similarity between the curves in figs. 7 and 12 in the 
thickness range 0 to 100 A. is indeed striking, the two series of measure- 
ments being in harmony. Holden’s suggestion that the reflection fringe 
asymmetry and the glass-silver reflection phase change are intimately 
related therefore receives experimental confirmation. 
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The experimental results may now be summarized and compared. 
with the predictions of Table IT. 

(1) For thick and very thin silver films the phase term y is close to an 
odd multiple of 7, the light absorption is small and the ratio R/R’ nearly 
unity. In an interferometer constructed of such films the reflection 
minima and transmission maxima are almost in perfect coincidence. 

(2) In the immediate neighbourhood of the minimum in the glass-silver 
reflectivity curve the absorption and the ratio R/R’ are large and 
x deviates from (2N+1)z. The reflection fringes are then usually 
asymmetrical. In the blue, where the maximum absorption is relatively 
small, the deviation is less pronounced than in the green and yellow for 
which wavelengths x closely approaches both 2Nz and (2N-+2)m thereby 
producing nearly symmetrical reflection fringes, the minima of which 
coincide with the transmission maxima. 

(3) In this same transition range y should pass through the value 
(2N-+-1)z in both the green and the yellow. As the ratio R/R’ will then 
be very high, the reflection and transmission minima might coincide in 
an air-gap interferometer using such a film for the upper plate XX’ ; 
this is particularly likely when the reflectivity of the rear surface YY’ is 
high. 

The optical properties of sputtered and evaporated silver films have been 
shown to possess similar characteristics. Gold films, whether they are 
sputtered, Rouard (1937) and Goos (1937), or evaporated, Sennett and 
Scott (1950), display the same general features as those of silver, and it is 
probable that these films will also produce reflection fringe asymmetry. 
Furthermore, the writer has observed asymmetrical reflection fringes 
using evaporated aluminium films, which Sennett and Scott have also 
shown to possess a peak absorption at a certain film thickness. The 
conclusions reached in this paper are therefore possibly typical of 
a number of metallic films. 

Metal films have been regarded by Rouard and others as homogeneous 
layers of definable thickness and the principles of optical interference 
have therefore been applied to them. From the oscillating nature of 
the optical properties of the films it is then concluded that the optical 
constants vary with film thickness. A more satisfactory approach 
appears to be to adopt the ideas of Maxwell Garnett (1904, 1905), who 
regards the film as consisting of many regularly scattered globules, the 
optical constants of which are the same as those of the bulk metal. The 
sphere diameters, which are small relative to the wavelength of light, 
are not necessarily uniform, the characteristic of each film being the 
volume of metal spheres in unit volume of the film. Krautkramer and 
Sennett and Scott have demonstrated that this theory provides a 
qualitative behaviour for many of the properties of metallic films, but 
that a number of features still remain unexplained. The theory might be 
improved by considering the non-spherical nature of the metallic 


particles. 
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By R. Firtu, 
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[Received September 1, 1950.] 


ABSTRACT. 


The notion of the velocity of the random movement of a particle and 
the problem of its observation is discussed in view of certain incorrect 
statements in a recent paper by D. K. C. MacDonald. 


THE problem of how to define and to measure the velocity of the 
random movement of a particle has frequently been discussed in the past. 
Recently D. K. C. MacDonald (1950) has contributed to the problem by 
applying to it the modern notion of the “‘ frequency spectrum” of a 
random movement and by deriving an expression for the ‘“‘ observable ” 
mean square fluctuation 82 of the velocity which he attempts to test on 
some early observations of the Brownian movement of small particles 
by F. Exner (1900). However, his analysis rather confuses the real issue 
and the application to the observations of Exner is partly erroneous. 
I, therefore, wish to make the following brief remarks on the problem. 
(1) Let x be a “unidimensional stochastic function” of time ¢. This 
means that a function x(t) does not exist in the usual mathematical 
sense, but that there exists a probability density function p(z,t) for 
observing a certain value x at time ¢t. Consequently all moments 2° are 
ordinary functions of time and their derivatives exist. In particular, 
V=d(x)/dt is the so-called drift velocity (which may be zero), and the 
quantity D=+ (x2)/dt, which, in general, is different from zero, is a 
convenient measure for the rate of fluctuations of x (D can, of course, 
be independent of time). But the quantity da/dt does not exist at all. 
If in particular x is a coordinate, a corresponding velocity component 
in the strict sense of the word cannot be defined. Vice versa, if a velocity 
component v is primarily introduced as a stochastic variable the 
corresponding position co-ordinate has no meaning. This is the well 
known situation in quantum mechanics where either coordinates or 
momenta of particles can be chosen as representing variables, but not 
both at the same time. It leads directly to the uncertainty principle 
for the simultaneous measurement of position and velocity. 
(2) There is, of course, no justification for regarding a physical 
quantity x as a stochastic variable under conditions where quantum 
effects can be neglected. For the random character of the variation of 
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x in such a case is only due to the influence of ‘‘ random external 
disturbances” on the observed system which are supposed to be not 
accessible to direct measurement. Thus x can be regarded as an ordinary 
function of time for an individual observed system, and the system, 
therefore, possesses at each moment a well defined velocity which, 
in the absence of a drift velocity, will irregularly fluctuate in time. 

One can, evidently, also define a “velocity” v* by the difference 
quotient 


over a finite time interval dt. This v* will be found to be practically 
identical with v as long as At is smaller than a certain time interval 7 
which depends on the mechanism of the observed system and which may 
appropriately be called its “ characteristic time ’’(T). 

If At<r, the statistical characteristics of v* must be identical with 
those of v. If, on the other hand, At is considerably larger than 7 the~ 
statistical characteristics of v* may be very different from those of »v. 
In any case, however, a statistical analysis of the behaviour of v* is 
absolutely equivalent to a statistical analysis of a record of displacements 
Ax at equal time intervals At, and so no useful purpose is served by the 
introduction of the quantity v*. 

(3) So far we have assumed that x is measured by means of an “ ideal ” 
instrument. In actual fact, the result of the observation by a real 
measuring device will usually be in the form of a graph representing the 
dependence of a “reading” y on time, y being a suitable measure of x. 
The rate of change w of y is obtained at any point of this graph by 
constructing the tangent to the arc of the curve at that point. This, 
of course, is only psssible with a finite accuracy, depending on the 
“resolving power ” dy of the instrument. 

The problem of the accuracy of the simultaneous determination of 
y and w has been discussed in detail by the author (1935, 1948), and he 
has shown that it is limited by a relation which is analogous to the 
quantum mechanical uncertainty principle. 


Within the limit dy of accuracy, one can also measure a quantity 
w* defined by 


wt 


A 

2 SY iy tie een eae. 
At ’ 

where At is again a finite time interval. Asa rule, it will be found that w* 
is identical with w as long as At is smaller than a certain time interval 
dt which we may call the “resolving power in time” of the instrument 
ee ee 


(t) 7 need not necessarily be the same as the so-called “relaxation time ”’ 
which plays a decisive part in MacDonald’s considerations ; this latter quantity 
is only defined for phenomena which are the outcome of the opposing influences 
of inertial and dissipating processes. 
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as the record can clearly give no indication whatsoever about the 
variation of the measured quantity x within time intervals smaller 
than déi(t). It has to be borne in mind, however, that both 8¢ and 
dy cannot be defined precisely and that their values are necessarily 
arbitrary within considerable margins. 

If one can make 6¢ of the order of magnitude + (or smaller) the record 
can be used for investigating the statistical properties of v. But if, as 
usual, d¢ is much larger than + the record can only be used for the 
investigation of the statistical properties of da for time intervals At 
larger than dt. This holds a fortior: for discontinuous observations of y at 
equal time intervals 4t. Clearly no complicated mathematical analysis 
is required to come to this almost trivial conclusion, and no analysis 
can reveal anything beyond the above statement. 

(4) Let us finally consider the practical example used by MacDonald 
to test his theoretical relations. The observations of Exner consisted 
in following, by means of a camera lucida, the movemnt of a Brownian 
particle as closely as possible with a stylus on a blackened glass plate. 
This curve was subsequently projected on a screen and the total length 
of the curve measured. The “ velocity’ of the Brownian movement 
was then taken to be this length divided by the total time taken for the 
observation. 

Evidently in this case d¢ is the “ reaction time ”’ of the observer for this 
sort of observation which is large compared with 7. Thus, according 


to (2), the measured “ velocity ”’ is simply (5s) /St where (5s) is the average 
are of the curve described in time 6¢. 

As the observed Brownian motion is a two-dimensional projection 
of the actual three-dimensional path, one finds easily for the probability 


density p(és,t) of ds in time t 


ry 2 
p(68.1)= sr exp(— Sr). tea eee 3) 


where D is the diffusion coefficient of a suspension of identical particles. 
* From (3) one obtzins at once 


6s)= J, ps, 9Bad08s)= Ve), Deir 


Thus finally 
is 041 ee eat () 


MacDonald, however, identifies his \/8? with w* which is incorrect, 
as his 8? refers to a one-dimensional Brownian movement; he also 


(+) St may have the character of a relaxation time but not necessarily so ; 
e.g. if the record is obtained by trying to follow the movements of a spot on 
a projection screen with a pencil, d¢ will mainly be determined by the time the 
nervous impulses take to travel from the eye to the brain and back to the 
muscles of the hand. 
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mixes up a mean absolute value with a root mean square value, a 
procedure he himself strongly criticises in his paper. 
For spherical particles of radius a 
fall 
— (ae . sib Ars . . . . . . . (6) 
Ta 


From (5) and (6) an expression for w* is obtained which, of course, is 
of the same order of magnitude as MacDonald’s expression (16), but there 
is no justification for his different numerical factor, although it is allegedly 
obtained by a refined mathematical analysis. 

From the values of w* given in Exner’s paper for particles of various 
sizes, one obtains from (5) and (6) values of dt between 0-1 and 0-2 sec. 
which are of the order of magnitude to be expected. No other 
conclusions can be drawn from these experiments. 
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SUMMARY. 


The “level” of integrated atmospherics received on narrow-band 
receivers tuned to a series of different frequencies between 3-5 and 50 ke./s. 
has been recorded. The band-width of the receivers was about 200 c./s. 
and the time of integration about one minute. The results were analysed 
so as to eliminate the unknown spectral distribution of energy at the source, 
and so as to show up the way in which propagation conditions changed 
with distance and with time of day. The results are shown in fig. 4, 
from which the most significant effect is seen to be a marked attenuation 
of waves of frequency below 10 ke./s. when travelling over distances of 
1000 km. or more. This attenuation is discussed and is related to the 
quasi-sinusoidal waveform often assumed by atmospherics which have 
travelled great distances. 

A decrease in the level of atmospherics has been found to occur during 
sudden ionosphere disturbances (8.I.D.’s) on frequencies below about 
10 ke./s. This contrasts with the known increase on frequencies above 


about 20 ke./s. 


§1. INTRODUCTION. 


A stupy of the propagation of the very longest radio waves entails the 
use of powerful sending stations, using large zrials, since the difficulty of 
radiating appreciable power increases as the wavelength is increased. 
For this reason the senders which have been used are those available for 
commercial radio services, and the longest waves at present available in 
this way have a length of about 18 km. (16 ke./s.). Much useful infor- 
mation has been gained about the ionosphere by studying the waves 
emitted by these stations (Best, Ratcliffe and Wilkes 1936, Budden, 
Ratcliffe and Wilkes 1939) and it is desirable that observations should be 
made on even longer wavelengths. The construction of special experi- 
mental senders for this purpose would be a large and expensive undertaking 
and an attempt has therefore been made to use instead naturally occurring 
atmospherics as sources of very long waves. This paper gives some of the 
results obtained at Cambridge in this way during the years 1948 and 1949. 

It is generally agreed that radio atmospherics have their origin in the 
electrical discharges associated with thunderclouds. The discharges, 
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being impulsive in nature, radiate over a wide frequency range and can be 
used as senders for investigating the propagation of very long waves. The 
large energy available in atmospherics at these wavelengths renders them 
easily observable compared with man-made and circuit noise. One of the 
major difficulties encountered in the use of atmospherics in this way is that 
the “ waveform ” of the original discharge is not known, so that the power 
radiated on different frequencies by any one atmospheric is not known, 
and the experiments must be planned so that the results are either inde- 
pendent of, or give us information about, the distribution of radiated 
power amongst the different frequencies. 

Another difficulty is that atmospherics are produced at different 
distances from the receiving point, and without some knowledge of, or 
assumptions about, their location, it is impossible to make useful deductions 
about the propagation of the associated waves. Fortunately we have been 
able, thanks to the assistance of the Director of the Meteorological Office, 
to make use of the chain of widely-spaced radio direction-finders which the 
Meteorological Office use to locate lightning flashes in their so-called 
‘ Sferics ’ organization. In this way it was known when fairly localized 
storm-centres were in action at distances less than about 1500 km. 

The method employed was to determine four points on the “ spectrum ” 
of the integrated received atmospherics, by recording the level on four 
different frequencies, and then to investigate how this spectrum changed 
from day to night, and with the distance of the storm centre. The 
assumption was then made that changes in the spectrum corresponded to 
changes in the propagation conditions, brought about either by the change 
from day to night, or by the change in the distance of the storm. 


§2. MerHop oF OBSERVATION. 


The results to be described were obtained from records of the average 
field strength of atmospherics on a number of different frequencies. The 
equipment used on one frequency is shown schematically in fig. 1 (a). 
A vertical «rial was connected to the common input of the various tuned 
receivers. Hach receiver selected the appropriate frequency component. 
from the atmospherics arriving. Its output was then rectified and 
averaged with a time constant of about one minute. The resulting 
voltage was recorded on paper moving at one inch per hour. 

The receivers were all of similar type. Each consisted of a tuned input 
stage followed by two tuned amplifier stages, and had an overall bandwidth 
between half-power points of approximately 200 ¢./s. The radio frequency 
output fed the detector circuit and was also available for viewing on an 
oscilloscope. The detector circuit is shown schematically in fig. 1 (b). 
The time constant for discharge of the condenser C was about 80 seconds 
which was very long compared with the average interval between the 
arrival of atmospherics. Hence a fairly steady output was obtained which 
represented a running mean of the input signal. The voltage Vp across 
the condenser C at any instant therefore depended on the signals received 
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during a time of the order of 80 seconds. In the simple circuit of fig. 1 (b) 
this voltage would affect the bias on the detector diode. To prevent this 
a voltage V, equal to Vz was connected in series with the diode. This was 
actually achieved, as shown in fig. 1 (¢), by applying the voltage V, to a 
cathode follower valve, and using the output of this valve to supply bias to 
the diode. The cathode voltage of a cathode follower exceeds the grid 
voltage by a small amount V, which is roughly constant, so that the output 
voltage of the cathode follower was Vg+Vx. The excess voltage Vx 
was compensated as explained below. 

When the radio frequency output was examined on an oscilloscope, it 
was found to consist of atmospherics of short duration and fairly large 
amplitude and of continuous “ noise ” of small amplitude, much of which 
was known to be due to ‘‘ man-made ” interference. On some frequencies, 
also, a small output due to signals from radio stations was sometimes 
observed. To reduce the effect of these unwanted voltages, an additional 
small voltage, V, in fig. 1(b), was connected in series with the diode, 
which prevented it from conducting until the level of the signal output 
from the receiver exceeded a small threshold about equal to V,. The 
voltage V, was actually derived, as shown in fig. 1 (c), from a constant 
voltage source V, which served also to compensate for the voltage Vx 
in the output of the cathode follower, so that V.=V)—Vx. Vz». was equal 
to 2 or 3 volts, which greatly reduced the contributions from the unwanted 
signals without appreciably affecting the wanted atmospherics. 

The apparatus was set up at the Radio Field Laboratory in Cambridge, 
where the general level of man-made noise was known to be low. Records 
were made between February 1948 and December 1949, at first on two 
frequencies, and after June 1949 on four frequencies. The frequencies 
were selected from time to time from amongst 3-5, 5, 7-5, 10, 14, 22, 27, 33, 
40 and 48 ke./s. The receivers were calibrated from time to time, at first 
by means of a C.W. signal generator, and later by means of a special 
impulse generator which produced a regular succession of pulses each of 
duration 1, 4 or 14 microseconds. These calibrations were only used 
to check the linearity of the receivers, and in the deductions to be 
described no use is made of the relative calibrations on the different 
frequencies. 


§3. RESULTS OBTAINED. 


The records were analysed so as to bring out effects which might be 
produced by changes in the conditions of propagation. For this purpose 
it was assumed, as a first approximation, that waves of all frequencies were 
propagated without appreciable attenuation from storm centres which 
were comparatively near (less than about 600 km.). It was next assumed 
that, in the absence of near storms, the atmospherics represented sources 
at a considerable distance (greater than 1000 km.) and that the recorded 
integrated levels measured relative to their values at the source would 
vary with frequency in a way which depended both on the nature of the 
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source and on the propagation characteristics. The effect of the propa- 
gation characteristic at night was established by comparing the normal 
night-time level with the level produced by near storms, as explained in 
sub-section (b) below. If, next, propagation conditions change as we pass 
from night to dayit should be possible to establish the change by comparing, 
in the absence of near storms, the levels recorded by day and by night at 
the different frequencies. This comparison is made in sub-section (a) 
below. It is convenient to consider these two comparisons in the reverse 
order. 


Fig. 2. (a-b) 
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(a) Typical atmospherics record in the absence of near storms. 
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(a) Change from day to night. 


Fig. 2 (a) shows a typical record for a 24-hour period when it was known 
that there was no storm centre within 1000 km. of Cambridge. From a 
record of this type it is possible to estimate, with some significance, the 
average level in the few hours before sunrise and the level in the few 
hours after sunrise, and to determine their ratio. The average value of 
this ratio, deduced from a large number of records obtained on different 
frequencies, is exhibited as a function of frequency in fig. 2 (0). Often 
jocal storm centres would develop during the afternoon so as to complicate 
the records, but by determining the daytime level from the few hours just 
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after sunrise, and at the same time disregarding the spasmodic increases 
produced by isolated near discharges, a significant ratio was obtainable. 
Fig. 2(6) appears to indicate that the onset of daylight conditions 
represents a marked increase of attenuation for frequencies less than 
10 ke./s. and greater than 20 ke./s. but near 10 ke./s. there is little extra 


attenuation. 


Fig. 3. (a-b) 
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(a) Atmospherics record, showing effect of near storm. 
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(b) Comparison of near and distant storms. 

Additional evidence is derivable from the increase in the average field 
strength which occurs when a comparatively close storm centre becomes. 
active. The type of record then obtained is shown in fig. 3(a). On 
several occasions records of this type were obtained, and it could be seen 
from the “ Sferics ’ charts that a comparatively local storm centre had 
become active at some distance less than 1000 km. Under these condi- 
tions it was possible to estimate the increase in average field strength on 
the frequencies being recorded. The results were studied by computing 
the ratio (average night field strength)—~(increase due to storm). The 
absolute value of the ratio depends, of course, on the intensity of the parti- 
cular storm being considered ; what is of interest, however, is not the 
absolute magnitude of the ratio but the way it changes with frequency. 
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For any one storm we have therefore normalized the ratios obtained on the 
different frequencies so that the ratio for a frequency of 22 ke. /s. is put 
equal to unity. During the short period of recording on 3-5 ke./s., only 
three storms of short duration and completely localized within 600 km. 
were available, and these provided the mean “ night-to-storm ”’ ratios for 
3-5, 5 and 22 kce./s. shown in fig. 3(b). Other storms, more extensive 
in range or duration, yielded values showing the same general trend at the 
low frequency end, while records made at different times on other fre- 
quencies yielded the points shown on the curve for frequencies between 
10 and 50 ke./s. In interpreting these results it should be remembered 
that the field on 3-5 ke./s. is predominantly the radiation field at distances 
greater than about 20 km. and all the storms considered were further away 
than that. 


§4. Discussion oF RESULTS. 


We start our discussion of the results by considering the implications of 
fig. 3(b). We first make the assumption, previously mentioned, that all 
frequencies are propagated equally well from a near storm centre, and that. 
the normal night-time level represents storms at a great distance. Then 
fig. 3(b) gives a measure of the relative attenuation suffered by waves of 
different frequency as they travel over great distances at night. The 
results are normalized so that the attenuation factor for a frequency of 
22 ke./s. is taken as unity. Of course the distances of the sources at night 
are very indefinite, but the results appear to show that, when the sources 
are at a considerable distance (greater than 1000 km.), the lower fre- 
quencies are much more heavily attenuated than the higher and that the 
‘“‘ low-frequency cut-off’? occurs just below 10 ke./s. Although this 
conclusion rests mainly on observatious made on a few occasions only, we 
consider it to be quite sound because it is supported by observations made 
on several other occasions when different frequencies were in use, or when. 
the storm centre was not accurately located. 

We next consider the curve of fig. 2 (b) which shows the extra attenuation 
introduced on the different frequencies as we pass from night-time to 
daytime. Since this curve refers to conditions when there were no storms 
within 1000 km. we assume tnat it represents propagation over great 
distances. Fig. 2 (b) is the result of a preliminary analysis covering the 
midsummer months of 1949 (more than 60 days’ observations are involved 
in the determination of the value for 5 kc:/s.), and we consider that it 
represents a well-established result and shows that the incidence of 
daylight has the effect of producing a considerable increase in the attenu- 
ation of waves propagated over great distances on frequencies below 
about 10 ke./s. | 

Fig. 4 (a) shows the attenuation curve of fig. 3 (b) for night, re-plotted 
on a logarithmic basis. Fig. 4 (b) represents the product of the curves of 
figs. 2(b) and 3(b), plotted in a similar way, and therefore represents the 
attenuation curve for day. 
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Although there may be some doubt about the absolute magnitudes of 
the attenuation factors shown in this figure, we feel reasonably confident 
that it indicates, in a general way, how the attenuation factor depends on 
frequency for transmission over distances of 1000 km. or more. We shall 
here concern ourselves with the results represented by the portion of the 
curve below 10 ke./s. 

In considering these results it is convenient to think of the space 
between the ground and the ionosphere as forming a four-terminal network, 
with the input terminals at the lightning flash and the output at the 
receiving erial. Our experimental results provide a description of the 
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Attenuation factor as a function of frequency for (a) night, (b) day, (c) day during 
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overall frequency characteristics of this network and, in particular, show 
that there is a sharp cut-off at frequencies below about 10 ke./s., but do 
not indicate how these characteristics come about. It might, for example 

result from the combination of separate ground- and sky-waves. In this 
connection, however, one point of significance emerges. The attenuation 
of a ground wave of frequency 5 ke./s. by absorption and diffraction, when 
travelling over a distance of 2000 km., can be calculated and the attenu- 
ation factor is found to be about 0-3. But the daytime attenuation factor 
measured on this frequency is much lower, about 0-02, as shown in fig. 4 (b). 
We deduce that the strength of the daytime field is much less than that of 
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the ground-wave alone. In explaining this result it will probably be 
necessary to remember that the wavelength (60 km.) at this frequency is 
comparable with the distance between the ground and the “ reflecting ” 
region, so that it is doubtful whether we are entitled to think of the 
problem in terms of separate ground- and sky-waves. 

It is not the purpose of this paper to discuss the explanation of the results. 
here presented, but I understand that Dr. K. G. Budden is preparing a 
paper in which he suggests a theoretical explanation. 

Before leaving the discussion of the curves of fig. 4 it is of interest to. 
relate them to the waveform of an atmospheric received on an aperiodic 
amplifier. If we assume that the atmospheric is a sharp impulse lasting 
considerably less than 210-4 sec. (the time period corresponding to a. 
frequency of 5 kc./s.), then simple transient analysis tells us what output 
we should expect if it is applied to a four-terminal network having the 
transmission characteristic of fig. 4. The result would be a quasi-sinusoidal 
oscillation with ‘* period ” roughly equal to the cut-off frequency of the 
network. 


Fig. 5. 
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Waveform of atmospheric from a distant source. 


Now it is well-known (Appleton and Chapman 1937) that atmospherics 
received from the greater distances very frequently show a smooth 
quasi-sinusoidal waveform with a period of the order of 150 microseconds. 
An example of one of these waveforms is shown in fig. 5 and it is clear that 
since it is so smooth and regular it would be difficult to interpret it in terms 
of a series of separate reflections from the ionosphere. It is, however, very 
much what we should expect as the result of applying an impulse to a 
transmission system with an attenuation curve as shown in fig. 4, and it 
seems probable that, once this attenuation curve can be explained, the. 
origin of the quasi-sinusoidal atmospheric waveforms will be clear. 


§5. Errecr or SuppEN JonospHERE DisTuRBANCES (S.I.D.’s). 


Bureau (1937) has shown that the Sudden Ionospheric Disturbances: 
(S.1.D.’s) of the type associated with solar flares, and frequently with 
‘“ Dellinger Fadeouts”, very often coincide with an increase of the average 
field strength of atmospherics recorded on frequencies in the range 
20-40 ke./s. We have confirmed his observations in this range, but our 
records on the lowest frequencies have revealed the unexpected pheno- 
menon that there the average field strength often decreases during the 
SI.D. This effect has been observed most frequently on 5 ke./s.. 
Examples are shown in fig. 6. 
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We have already seen that, during the day, the average field strength 
on the lowest frequencies is very small, unless there is a fairly near storm, 
and indeed if there are no storm centres within 2000 km. the average field 
strength recorded on 5 kc./s. is so low that a further decrease due to an 
S.I.D., if it does occur, is not noticeable. For a frequency of 3-5 ke./s. the 
corresponding distance is about 1400 km. 

If the storm centre is at distances between 1500 and 2000 km., it has 
been observed that S.I.D.’s can reduce the mean field strength by a factor 
of 10:1 on a frequency of 5 ke./s., and for distances between 600 and 
1500 km. a reduction of 3:1 has been observed. Observations on 
3:5 ke./s. have been less extensive, but on one occasion the field due to a 
storm centre at a distance of 800-1500 km. was reduced to one fifth during 


an S.L.D. 
Fig. 6. 
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@ 8 October 1949 


Records of integrated atmospherics at times of Sudden Ionospheric Disturbances. 
(a) 8th October, 1949 (zero level shown only for 3-5 ke./s.); 
(6) 6th November, 1949. 


Although the effects associated with S.I.D.’s have depended on the 
distance of the storm centres, we have attempted to represent them in a 
general way by the curve (c) in fig. 4. In examining this curve it is 
interesting to notice the different effect of an S.I.D. on the resultant signal 
produced at a great distance on different frequencies. For frequencies 
below 10 kc./s. the signal decreases : for frequencies between 10 ke./s. and 
about 200 ke./s. the signal increases : for frequencies above about 1 Mc./s. 
the signal decreases (Dellinges effect). A satisfactory theory of propa- 
gation will have to explain these results. 
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SUMMARY. 


When a body travels through a gas containing an isotropic distribution 
of sound waves it suffers a retarding force proportional to its area, to: 
the energy density of the sound, and to its velocity. Three effects might 
contribute to this force: (a) the distribution of sound with respect to 
the moving body is no longer isotropic ; (b) the apparent wavelength of 
the sound is altered by the motion of the body ; (c) the non-linear terms. 
in the equation of motion lead to a coupling between the sound waves 
and the steady flow of gas round the body. Effect (c) is shown to produce 
no force. Effects (a) and (b) are evaluated in the case of a sphere, and 
the force is shown to be expressible in terms of the force exerted by a 
single sound wave on a stationary sphere, already calculated by L. V. King. 


$1. INTRODUCTION. 


Iv is well known that a body moving uniformly through a non-viscous 
fluid at rest is not acted on by any force. The resultant couple also. 
vanishes if the body is symmetrical (Lamb 1932). If the body is at rest, 
and a plane sound wave impinges on it, the sound wave is scattered, and 
a force is exerted on the body. This force has been calculated by Rayleigh 
(1905) in the case where the obstacle is the flat end of a cylinder containing 
the gas, and by King (1934) for the case of a spherical obstacle. (King 
actually treats the more general problem of a free sphere of density p. 
The result for a fixed sphere is obtained by taking p infinite.) An 
apparent discrepancy between their results is discussed in §2 of the 
present paper. 

If a symmetrical body hes at rest in a gas containing an isotropic 
distribution of sound waves, the resultant force on it clearly vanishes. 
If it moves with velocity v, it experiences a retarding force proportional 
to its area, to the energy density of the sound waves, and to v/c, where 
c is the speed of sound. In §3 a plausible argument is advanced for 
believing that this force should vanish. The exposure of this fallacy 
illustrates an unexpected difficulty which arises if the thermal motion 
of the molecules of a gas is analysed into Fourier components in the way 
customarily used for solids. The result to be obtained is proportional 
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both to the speed of the body and to the energy density of the sound 
waves, and therefore involves terms of the third order of small quantities 
in the equations of motion. In § 4 the equations are developed to the 
necessary order. In §5 the explicit interaction between the uniform 
motion of the body and a single sound wave is evaluated and found to 
lead to no resultant force. However, an isotropic distribution of sound 
waves in the gas is not isotropic in axes attached to the moving body, 
and in § 6 it is shown that this anistropy leads to a force of the kind 
expected. 

This problem has a certain intrinsic interest. I was actually led to 
study it by reading a paper by Leibfried (1950) in which he applies similar 
considerations to estimate the resistance to the motion of a dislocation 
in a crystal caused by its scattering of sound waves in that solid. His 
result for this more complicated problem depends upon several drastic 
simplifying assumptions. He finds a resistance so large that the terminal 
velocity of a dislocation in a crystal under a reasonable stress could never 
become comparable with that of sound. Since this result would render 
impossible a type of process which plays an important réle in some modern 
theories of plastic deformation (Frank 1948), it seemed worth while to 
treat this related, but simpler, problem in detail. The results are 
compatible with that of Leibfried. 


§ 2. Discussion oF RAYLEIGH’sS AND Kine’s FORMULA FOR THE 
PRESSURE OF A SOUND WAVE. 
King considers a single plane wave of wavelength A incident on a 
sphere of radius a. If the energy density of the sound wave is E, the 
force P on the sphere is given by 


eV) eee Oe fee ae eee (1) 
where i ANC es gee BY gee a edt ) 


and it is easy to see that the function L must tend to the value 3 When. 
ka is large. King’s computation leads to L(1)=0-425. The distribution 
of scattered radiation from a large sphere may be resolved (Morse 1948) 
into a forward lobe which interferes with the incident beam to form the 
shadow and an isotropic reflected beam. Thus, of the radiation inter- 
cepted by the scattering cross-section 27a”, half continues in the forward 
direction with unchanged momentum, and half is scattered isotropically 
with zero total momentum. If the momentum in unit volume of the 
incident wave is m, the change of momentum in unit time caused by the 
sphere is mcza*, and comparison with (1) leads to 


elle EP cos 6 Sea CD] 
This appears to disagree with Rayleigh’s result. Both authors assume 
that the pressure p in the gas depends only on its density p according to 


the law 
UE, ee 
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If the density of the undisturbed gas is po, the velocity of sound c is 
given by 


ctefi(p,)it 6 et OT ee re (5) 
Rayleigh’s formula for the momentum is, in the present notation, 
m==(L-t-pgf (pp) 207/26. aeestieeenss mae ee Bale) 


The discrepancy between (3) and (6) does not arise from an implied 
assumption f”(p))=0 (corresponding to Boyle’s Law) in the derivation 
of (3), but from a difference in the physical situations considered. A 
similar ambiguity is known (Andrejew 1940) to arise in the determination 
of the energy density in a sound wave. King treats a finite obstacle 
in an infinite gas. Because the plane sound wave is disturbed by the 
obstacle, the time average of the density of the gas varies from place to 
place, and this variation contributes to the resultant force. Rayleigh 
considers a fixed quantity of air in a cylindrical pipe, either at rest or 
containing a resonant sound wave in a quasi-stationary state. The pipe 
is closed by a piston at one end and by a rigid wall at the other. Rayleigh’s 
argument will be developed here in the notation which will be convenient 
in later work. 

In Cartesian coordinates x,;, the components of velocity u; are given by 
the equations of motion : 


Du; op 
p Di Oe Cie? : : . . . . . . : (7) 
Deare 0 
where Dieronn Mt Oa,” A . - < 4 : = ‘ (8) 


The repeated suffix 7 implies a summation over the values i=1, 2, 3. 
If a function a is introduced by the relation 


dp 
eas Ree! ko Sd ye Fameing 
the equations of motion become 
Du, Ow 
pte x (10) 
If the motion is irrotational, there is a velocity potential ¢ satisfying 
U=— _ (11) 


(Rayleigh here uses —¢.) 
The equations of motion are satisfied if 


_ 06 1/8d\2_ Dé  1/ad\2 
=F —3(5e) — pet 5(g) - ir RS ann bo) 


Changing only the notation, Rayleigh’s argument now proceeds: “If 
we integrate over a long period of time, 4 disappears, and we see that 
Jo dt+-}J(0¢/0x)? dt retains a constant value at all points of the cylinder. 
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The value at the piston is accordingly the same as the mean value taken 
over the length of the cylinder.’ The pressure on the piston is derived 
from the value of w at that point. According to (9), the definition of a 
involves an arbitrary constant of integration. Rayleigh determines this 
constant from his ‘“‘mean value” result just quoted, using his assumption 
that the total mass of gas, and hence the mean density taken over the 
length of the cylinder, is unaltered by the presence of the sound wave, 
The term in Rayleigh’s formula for the pressure on the piston which 
depends only on conditions near the piston would lead to (3). The term 
in f” in (6) arises from integrals over the length of the tube, expressing 
the constancy of the total mass of gas. 

Two points in the argument quoted deserve rather more attention. 
Firstly, it is obvious that ¢ does not vanish in the time integral. It is 
necessary to take the time average before ¢ can disappear. If ¢ actually 
did ‘‘disappear’’ in the integration, the time average of w+ 4(d¢/dx)? 
would not be a determinable constant but zero, leading again to (3). In 
fact the arbitrary constant of integration in w corresponds to an arbitrary 
addition to ¢ of the form t/(x). This gives an aperiodic contribution 
—t dy/dx to u,, and the postulate of a steady state leads to u(~)=constant, 
which justifies Rayleigh’s method. 

The second point is that when terms of the second order in the wave 
equation are considered, a true stationary state almost certainly does not 
exist (Rayleigh 1910). Time averages can only be given a meaning by 
a special process. If the unperturbed motion is represented by a 
sufficiently small periodic function ¢5, the second order terms are small 
perturbations. After time t, ¢ becomes ¢,.+-5¢, where 6¢ is in general 
of order ¢2¢. If 5¢ in fact has this form, the time average of d¢/dt does 
not tend to zero, but tends to an unknown quantity of order 43, which 
is of the same order as the terms to be evaluated. A periodic state may 
be maintained near the piston by using a train of waves in an infinite 
tube, but this no longer requires that the mean density should be unchanged. 
Alternatively, the rigid end wall of the cylinder may be replaced by another 
piston which is controlled from outside in an attempt to maintain a 
periodic state inthe tube. This artifice can only be used if it can be shown 
to be possible, and in particular it must be shown that the piston will not 
continuously supply energy to the gas in the tube. The second order 
terms in the wave equation lead to an irreversible change, and thus 
entropy is continuously generated in the gas. In a steady state this 
entropy flows out through the walls of the tube as heat, while the 
controlled piston supplies energy at the same rate. This is exactly the 
experimental arrangement which would normally be adopted to measure 
the pressure of radiation in a closed tube. It is remarkable that a constant 
flow of mechanical work into the tube and of heat out is not a result of 
experimental imperfections, but is essential to aself-consistent interpretation 
of the results. 
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§3. RELATION BETWEEN THE RESULTS OF THIS PAPER AND THE KINETIC 
THEORY OF GASES. 


Equation (1) gives the force on a sphere at rest in a single sound wave. 
If there is an isotropic distribution of sound waves of the same wavelength 
and of energy density E, the force on a sphere at rest clearly vanishes. 
It will be shown in §6 that if the sphere moves with velocity v it 1S. 
subject to a retarding force P given by 

P=na*H(v/c)Li(ka), . . » » - = =. (18) 
where L(kay=4 (40 (ka)? kalikay| ee ea 
The first term in (14) arises from the ‘‘ aberration ”’ of the sound waves, 
the moving sphere meeting a larger flux of sound waves on the leading 
than on the trailing hemisphere. The second term arises from the 
Doppler effect, the waves incident from ahead having a shorter effective 
wavelength and a correspondingly larger scattering cross-section than. 
those incident from behind. 

Consider a large sphere moving through a monatomic gas at rest. 
The pressure of the gas is caused by the motion of the gas atoms. In 
the case of a solid it is customary to resolve the thermal motion of the 
atoms into normal modes of vibration of the lattice, and it is plausible 
to do the same in the case of a gas. In both cases it is necessary to consider 
vibrations of wavelengths as short as twice the interatomic spacing, and 
nevertheless the vibrations are regarded as those of a continuous medium. 
These short waves are far more numerous than the larger waves, for which 
such an approximation is reasonably accurate. In both cases the scattering 
of short waves is so large that their mean free path is of the order of one 
wavelength. Since this theory is reasonably successful for solids, it 
might be expected also to succeed for gases. If the whole thermal 
energy of the gas is resolved into sound waves in this way, the moving 

sphere is retarded according to (13), with kaS>1. If the pressure of the 
gas is p, then Pe while L(00)=%. So 


P=ra2poic. J) cider oy satires 


A retardation of this enormous value would leave billiards as the only 
practicable ball game. While the approximations leading to (15) could 
alter the value of P appreciably, it is obvious that the derivation is not 
merely inaccurate, but basically wrong. 

The same problem may be treated in a more orthodox way by regarding 
the pressure of the gas on the sphere as being produced, not by the 
impact of sound waves, but by the impact of atoms. Again, more atoms 
strike on the leading hemisphere than on the trailing, and, if the sphere is 
smooth, formula (15) is again obtained. In this case the error in the 
assumptions is obvious. » The result can only be expected to hold if a is 
much less than the mean free path of atoms in the gas. The distribution 
of momenta among atoms striking an element of the surface of the sphere 
is then the superposition of two distributions, corresponding to incident 
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and reflected streams of gas. If a is large compared with the mean free 
path, the distribution of momenta among atoms reaching the surface is 
not that corresponding to distant parts of the gas, but corresponds to 
the pressure and mean velocity of the gas (which is parallel to the surface) 
at that point in the ordinary hydrodynamical theory. Forces of the 
kind just considered arise only when v changes, and give rise to an increase 
in the effective mass of the sphere. 

In the original derivation of (15), where the thermal motion of the gas 
was resolved into sound waves, it is now clear that the short waves, 
which contribute most of the pressure p, do not contribute to P because 
of their short free path, but there remain enough long waves of very long 
free path to give a retardation far greater than is observed. The basic 
difficulty in regarding the pressure of a gas as the sum of the pressures 
of all the sound waves in the gas now appears. It follows from (9) and (12) 
that in any standing wave the time average of the pressure is greatest at 
the nodes, where (0¢/0x;)?=0, and least at the antinodes. Since the 
wall is a node for every wave, the pressure on the wall exceeds the pressure 
in the interior, which is equal to the ‘internal pressure” of the gas, 
that is, zero for a perfect gas. In a solid, this resolution of the thermal 
motion into independent sound waves is an acceptable approximation, 
because each sound wave can be propagated alone through the elastic 
solid. This is not truein the gas. The velocity of propagation of a sound 
wave is proportioanl to T!, where T is the absolute temperature. If the 
whole thermal energy of the gas is resolved into sound waves, T is a measure 
of the number of sound waves already present, and the propagation of one 
sound wave depends essentially on the presence of others. To analyse 
the motion into independent sound waves is therefore meaningless. If, 
in the problem of a sphere moving through a gas at rest, the thermal motion 
of the gas is resolved into sound waves, the disturbance of one particular 
sound wave by the sphere is produced not only directly, but also 
indirectly by the interaction of this sound wave with other sound waves 
already disturbed by the sphere. In the hydrodynamic solution, the 
points of stagnation are points of high pressure and so, in adiabatic flow, 
of high temperature and of high sound velocity, while the equator is a 
region of low sound velocity. This temperature distribution refracts 
sound waves so that they tend to travel along the lines of flow, in the 
same way that, in the molecular picture, incident molecules are deviated 
by collisions towards the lines of flow. Thus (13) cannot be disproved 
by deducing from it the absurd consequence (15). 


§ 4. THe Basic Equations. 


The motion of the gas is governed by the equation of continuity, 


O(pu;) _ 15 
= + ax, = 0, a") Vébe 8 0, eR Ley ( 5) 
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which may be written 
Dinp &¢ 
hee Be Jn i iawvetr nc Ieee? 


It follows from (4) and (9) that p may be expressed as a function of & 
alone, and then from (12) that p may be expressed in the form 


PEeO(O) 40 or ees a ae ee ne 


The function p(¢) involves both ¢ and its derivatives. The equation of 
motion (16) may now be written in terms of ¢ alone as 


Din p(g) _ 4 
em remeat rite ey a, * (18) 


while the pressure is given by 


=f pl) lie ve ae, tes cee EL 


These equations will now be expanded in powers of ¢, using as an 
auxiliary variable the condensation s, a pure number of the same order 
as ¢, defined by 


p= po\1--8). bose h ary ~ ee es chs, age ae 

To the third order (4) becomes 
P=Pot prcrs+4pef’se-+ apes, . 2 . . . (21) 
where pee ilps ie sekg ieee Oe lie Sonera: (22) 


and the derivatives f” and f’” are evaluated at p=p,. To the same order 
(9) becomes : 


w= C'8—4(O— pof")?+3(C— pof” +3p0f")s®. . . . (23) 
Inverting (23), 
ep Wee: 7 ow of” Up 2 wo 
$= BTR pol”) Te tal —4p00f" + 3pof"— poet”) sa. - (24) 


Using (12) and (24), (17) becomes 


Po 9 1 po (Od Lip my (Ob\? 
Pa ao + 5 (se) +5 5 (e—eaf") (FF) 


1 Pp A g my we 0 
J gto (4p 92h" +3p5,f"— paerf 1G; x) ee 546595 
It is convenient to introduce the abbreviations 


2 ion. 
ae = $5 i; i Ox. at = = $;; ae == Piz; etc. 5 (26) 
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With this notation the equation. of motion (18) becomes, to second order, 
© ;— $= —2$ 6 ;— por 4f"od, . . « . . (27) 

and (21) becomes, to third order, 
P=Pot Pob—bpobi+ bpoe°d?— 4 poe *bdi + bpoe (C>—pof")$%. (28) 
It is noteworthy that f’’’, which appears explicitly in (21), isabsent from (28). 


§5. A Bopy Movine wirn Unirorm Vexociry THROUGH A Gas 
CONTAINING A SINGLE PROGRESSIVE SouND WaAvE. 


In treating the interaction of a sound wave with a body moving with 
uniform velocity v, it is convenient to use axes (x{ x} 24) travelling with 
the body, and also to make a Lorentz transformation of the time. The 
new coordinates are given by 


(1—v?/c?)taj =x ,— vt, 
Le=Xo, 


(kami ICE), 


(il 
U3=%3, 


(1—v?/c?)*t’ =t— va, /C?. 

Ifv/cis treated as a small quantity of the first order, the wave equation (27) 
takes the same form, to second order, in primed coordinates. In (28) the 
term po? does not contribute to the time average of p, and the remaining 
terms are unaltered, to third order, by the change of coordinates. (By the 
argument of § 2, a term in ¢ proportional to ¢ can only lead to a change in 
pressure which is independent of position, and so produces no force on a 
closed surface.) Equations (27) and (28) may therefore be considered to. 
hold in moving coordinates. They are ‘“‘ mixed” equations, in which 
p and:¢ retain their original values measured in fixed coordinates, while 
‘the differentiations are carried out in moving coordinates. 

The present problem is concerned with the interaction between two 
disturbances. One, representing the uniform motion of the body with 
speed v, is independent of time ; the other is, in first order, a sinusoidal 
function of frequency w. The first order solution of (27) is then 
expressible in the form 


dV=V(ax;)+C(x,) cos wt+S8(a,) sin wt. 

In moving coordinates, the disturbance represented by C and 8 is still 
periodic, and the first order solution may be written 
$V=V[ (ai +t’ )(1—v?/c2)-4, x3, #5]+C'(a') cos w’t’+8'(w)) sinwl’. . (30) 
If the solution to second order is 4=4¢+4, ¢™ is a solution of the 
equation 
ge" — — $2" = —20'(V5+C' cos w’t’--8% sin w't’)(—C', sin w’t’ +S, cos w’t’) 

ron Pe '3(—C’ sin w’t’-L8’ cos w’t’)(C’ cos wt’ sin w’t’). 

(31) 
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Here ¢/, is an abbreviation for 0°4/dx;”, and C; is an abbreviation for 
aC'/dx’, etc. Terms in 0V/dt’ in (31) are of the third order, and have 
been neglected. 

All terms in 4”) are of frequency 2w, except those arising from the 
term 2w/V‘(C, sin w't/—S', cos w’t’), and the solution for ¢ may be written 
to second order as 


d=h044 
=V-+C’ cos w’t' +8’ sin w’t/ +C®” cos 2w’t/+S8@” sin 2w't’ 
+(VC) sin w't’+(VS) coswt’, ©... .. . Seine iia >)) 
where (VC)’ is a solution of 


VO) -fw?(VC) S20 iC, eas 


and (VS)’ is a solution of 
c2(VS);4 @2(VS)'=—20'VjS; . . -. 2 (84) 


When this solution is substituted in (28), and terms of higher order 
than the third are neglected, it is readily seen that the time average 
Pm of p is given by 

Pm=Po— tP0Vi— 4Po(C 7? +877) + dpe 2w'?(C?+8) 3 
+ 4polC'(VS),+8;(VC)/]+ Fp,e~%eo’2[C"(VS)-+8'(VO)'] 
FPG OV (OSS. Cy) Soe ae ee ae 
Here p, is the pressure of the undisturbed gas, the second term is the 
pressure due to the steady flow, which gives no force and, for a symmetrical 
body, no couple. The next two terms give the acoustic radiation pressure, 
as calculated by King, and the remaining three terms give the pressure 
produced by the interaction of the steady flow and the sound wave. 


Using equations (33) and (34), together with the fact that C and V 
satisfy equations of the form 


2C;,,+@2C=0, em MEA Wein eli) 
reduces the last three terms in (35) to the form 
—polC'(VSY +8'(VC)'], 2 2. ee. (87) 


It may be shown that this term gives no resultant force on the body. 
For a smooth function Q may be defined inside the body, which has on 
the surface of the body the same first and second derivatives as 
[C’(VS)'+S(VC’)]. If dd" is an element of area of the surface, the total 
force on the body in the direction a; is 


3p | [ Q d2j=—4po | | [ (2Qj:/22%) dr’, 
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where the integral dr’ is taken through the volume of the body. This 
may be written 


0 Fi ] eke 
—FPo an, {ff Q;,47’= — tpo Ox, iH Qa". 
At each point of the surface, 
Qi=C{VS)'+C'(VS),+8(VCy +8'(VC);, 
-and each of the derivatives in this expression is orthogonal to d2 , because 


-each term C’, (VS)’, S’, (VC)’ represents the velocity potential of a flow with 
the surface of the body as a fixed boundary. 


§6. A Bopy Movine THROUGH AN IsorrRopic Ftux or SounD WAVES. 


Consider a body moving with velocity v through a large space containing 
an isotropic flux of sound waves of frequency w. It follows from § 5 
‘that the third order terms in (35) do not produce a force on the body. 
However, the radiation is not isotropic in coordinates moving with the 
body, and, even in the case of a sphere or a cylinder, the pressures 
corresponding to the second order terms in (35) lead to a resultant force. 

In a fixed coordinate system a typical sound wave has the velocity 
potential 

d=dy COS (kX, +hot.+hytz—wt), . . . . . (38) 
where 
a se ee (rrr 22) 


Referred to axes moving with the body, as given by (29), the velocity 
‘potential is 


b=¢y cos [(k,—wv/c?)(1—v?/c2) “ar +- haa + hgx3—(w—kv)(1—v/c?)-7¢], 
oe aren) 


Neglecting terms of order v?/c?, this is a wave travelling in the direction 


Vee | ay a 
Ry «Ky kgv 
fs a De eee 25a (4D 
ls k + kc ? ? ( ) 
, hg , kykgv 
Lp k a kc ’ 
with frequency 
w'=w—kv (42) 
-and wave number 
esicr ah ta Ce We P43) 


Equation (i) may be written in the form 
P=4mp $r(ka)PL(ka), . . . » » «+ (44) 
SER. 7, VOL. 41, NO. 323.-—DEC. 1950 At 
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and by substituting (40) in (30) and (35) it is clear that the force on a 
moving sphere is given by (44) with k replaced by k’. The resultant 
force produced by the isotropic radiation must act in the direction 7,.. 
The force produced by the wave (k,k,k,) has a component in this direction 
of 


LamppPieayrlLak ay... 28 as) sun eLeOy 
On setting 1,=k,/k, this may be written, to first order in v/c, as 
drrp bh {l, k?a?L — (v/c)[(1+2)ka2L+earL’'}}. . . . (46) 


If ¢, d2 is now taken to represent the velocity potential of waves in 


a small solid angle dQ, the total force P on the sphere is obtained by 
integrating (46) over all directions as 


= 2 27.272 
pa 27 Potol erat (ka) + kal (ka). eee hia: 


The energy density in the sound flux is 
E=2rp,¢2k?, . . . . . . . . (48) 
and (47) may thus be written 


Ba TEP" [41 (ba)+ kal (ka)), 1 ee 


This agrees with (13) and (14). 
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CXIV. A Note on the Decay of Current in Germanium Diodes. 


By J. R. Trttman and H. Yeo, 
The Post Office Engineering Research Station, London*. 


[Received October 4, 1950.] 


A modulator using germanium crystal diodes, for use in time division 
multiplex systems, shows an effect which enables a property of the- 
crystal diode—not hitherto reported on—to be investigated. 

Fig. 1 (a) shows part of a multiplex system consisting of two channels 
only. (E,;),—A sin w,f is the source of modulation for channel 1. At 
all times other than when a PI pulse is applied the diode W1A passes 
about 1 ma. and has low impedance, R, (1502), to a.c.; the diode 
W1B is therefore subject to a “ reverse ” voltage of about 1 and presents 
a high impedance, R,, to a.c. When a P1 pulse is applied W1B passes 
a pulse of current whose value lies between 2-2 and 3-8 ma. according to 
the instantaneous value of A sin w,t; WLA becomes of high impedance. 
The modulator of channel 2 works in a similar way, with sin w,¢ for sin w,f, 
P2 for PL, etc. ; its output is connected in parallel with that of channel 1. 
The pulse length is fixed at 0-2 psec. but the relative timing of 
Pl and P2, ¢,, is controllable. The trains of pulses, designated HP1, 
HP2, are together transmitted to the receiving end of the system where 
they are separated and demodulated by low pass filters. The two outputs 
(E,), and (Ep), substantially reproduce the waveforms of (E;), and 
(H;), respectively, with equal efficiency. Each, however, contains 
components, at low levels, in addition to that wanted. These include 
a crosstalk component ; thus (E,), contains a component C,, in wet and 
(Ey). a component C,, in wf. 

The circuit is so designed that the crosstalk ratio C,,/E), expected, 
is small; the crystal diode W1B is chosen to have R,>1MQ in order 
that the component of C,, arising from the finite values of R, and R, is 
<E,/5000. Only for values of ¢,<0-5 psec. will the attenuation distortion 
at high frequencies, introduced by residual shunt capacitances in the: 
common path, cause C,, to exceed E,/5000. 

However, in practice, values of C,, of up to E,/100 have been noted. 
at t,=1-0 ysec., depending on the sample of diode used as W1B. These 
values decay approximately exponentially, falling to l/eth in about 
0-4 to 1-2 psec. (see fig. 2). The use of either a thermionic diode or some 
other samples of crystal diodes as W1B reduces the crosstalk to about the 
expected low level residual. 

The newly noted crosstalk voltage is always out of phase, by 7, with the 
signal obtained with ,=0, and hence with any leak due to the finite: 
value of R,. Any failure on the part of R, rapidly to reach its steady 
state value after the cessation of P1 is, therefore, not the cause of the new 
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crosstalk voltage. This conclusion is supported by experiments with 
the modulator of channel 1 changed to that of fig. 1(b). Here W1B is 
a thermionic diode, or a crystal giving no significant crosstalk when used 
as WIB in the circuit of fig. 1 (a). Crosstalk still arises when a crystal 
giving crosstalk in the position of W1B of fig. 1 (a) is used as W1D in 
fig. 1(b), despite the fact that signal leak is now adequately suppressed 
by W1A and WI1B. 


Fig. 1. 
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AOD - (a) Schematic diagram 
of the multiplex system 


(4) The inc Fed modulator 


The explanation of the new crosstalk seems to lie in the slowness of 
dispersal of the charge carriers. But the evidence available on transit 
times does not afford a simple explanation. 

Lark Horovitz (1946) has noted that the rise of current through 
a germanium diode is noticeable slower than the rise of an applied voltage, 
when the latter takes place in 10~’ sec. ; the decay of his voltage pulses 
was so slow as to obscure the effect now noted. Douglas and James 
(in press) suggest that the finite transit time of the charge carriers may be 
responsible for the decrease of rectification efficiency noted at frequencies 
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algae about 10 Mc/s. They have correlated this falling-off with increasing 

“ turnover ” voltage, V, and have pointed out that turnover voltage 
may well be a measure of the thickness of the barrier layer. They do 
not give a value to any transit time, but in view of the fact that the 
rectification efficiency is not grossly impaired at 50 Mce/js., a time of 
10-8 sec. rather than 10-® sec. seems to be indicated. 

Correlation between the crosstalk voltage and other properties of the 
crystal diodes was sought, even though the large differences of crosstalk 
voltage noted (see fig. 2) could be obtained from units differing little 
in V,. In general, the crosstalk voltage was smaller the lower the 
turnover voltage ; two batches of crystal rectifiers type CG5 supplied 
by Mr. Kinman of the B.T.H. Co., Rugby, having V, less than 80 volts, 
consistently gave C,,.<E,/2000 at t,=lysec., whereas some CG4’s 
(V;>100 volts) gave C,,—E,/200 at t,=I1yusec. At the same time some 


Fig. 2 
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The variation of C,, with time for four samples of diode. 


other CG4’s gave C,,<E,/3000. Units with low turnover voltages often 
have low values of R,. There seemed no correlation between C,, and 
R,, however ; in fact, one of the two batches of CG5’s giving small values 
of C,, had, in the main, values of R, of about 100k2, whereas the other 
batch was specially made to have values of R,>1 MQ. 

Some of these experiments point, therefore, to the effect being 
independent of the operation of the diodes which gives them their better 
known properties ; further work based on this indication is in hand. 
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Superconductivity of Niobium. 


By L. C. Jackson and H. Preston-THomas*. 
H. H. Wills Physical Laboratory, University of Bristol. 


[Received October 10, 1950.] 


Untit recently, published data on the properties of the “hard” 
«superconductors were both scanty and inconsistent, particularly in the 
ase of niobium. Studies have been made in this laboratory of the 
‘properties of the latter metal, both before and after heat treatment. 
._At the start of the investigation the only published measurements of the 
-critical magnetic field of niobium were those of Daunt and Mendelssohn 
(1937) in the liquid helium region. A cryostat was therefore constructed 
‘to extend the measurements to the region between 4:2° K. and the normal 
transition temperature, given by Meissner, Franz and Westerhoff (1933) 
as 9-2° KK. A description of this apparatus is to appear elsewhere. Three 
methods have been used to determine the superconducting properties ; 
‘resistance measurements, an A.C. inductance method identical in essentials 
with that of Webber, Reynolds and McGuire (1949) and the ballistic 
method of Keeley and Mendelssohn (1936). The satisfactory operation 
of the apparatus was checked by carrying out measurements on pure 
lead (H.S. Lab. no. 11,116) by the three methods mentioned above. 

The critical magnetic fields obtained by the ballistic method for a 
sample of niobium (Johnson, Matthey and Co. H.S. Lab. no. 1673— 
‘spectroscopic examination showed the presence of traces of molybdenum, 
nickel and silicon together with a very slight trace of aluminium— 
used as received, not heat treated) are given in fig. 1 ; curve A, beginning 
of penetration of magnetic field, curve B completion of penetration. The 
amount of frozen-in flux is not large, of the order 10 to 15 per cent and 
varies with the temperature. Curve C of fig. 1 is derived from the 
alternating magnetic field measurements and agrees closely with the 
‘curve determined from measurements of the electrical resistance of 
niobium tape rolled from the same sample. The spread of the transition 
‘curve in the absence of a magnetic field is about 0-07 deg. 
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Curves A’, B’ and C’ of fig. 1 give the corresponding results for another 
sample of niobium. (Johnson, Matthey and Co. H.S. Lab. no. 2088— 


trace of silicon and very slight trace of aluminium detected spectro- 


Fig. 1. 
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scopically—used as received, no heat treatment.) Ourve D, fig. 2 a 
the critical field curve for niobium recently determined by Cook, Zemansky 
and Boorse (1950) using the alternating magnetic field method. 
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Fig. 2 shows the results of the ballistic method for an older, less pure: 
specimen of niobium (Hilger H.S. Lab. no. 8290—spectroscopic analysis 
showed the presence of about 0-2 per cent tin and possibly slight traces. 
_of copper, iron and gold—used as received, not heat treated). This is 
also the Lab. no. of the sample used by Daunt and Mendelssohn (1937) 
in their determination of the critical magnetic field between 1-5 and 4-2° K. 
Their measurements showed that about 20 per cent of the critical field 
remained frozen in the specimen. The solid line A in fig. 2 gives their 
results. The dotted line B is the interpolation between 4-2° K, and the: 
normal transition temperature determined by Meissner et al. (1933) as 
given in Shoenberg’s ‘‘ Superconductivity ” (1938). The points on the 
dotted line C are those at which the penetration of the magnetic field 
into the specimen is complete. Alternating magnetic field measurements. 
with the same specimen gave a normal transition beginning at approxi- 
mately 8° K. and becoming complete at approximately 7° K. The 
spread of the transition appears to become smaller as the applied field 
is increased but it is still too great to provide useful data for the H,—T 
curve. The material is, moreover, so brittle that it is impossible to roll 
tape specimens from it suitable for resistance measurements. 

There is a considerable similarity between the above results for niobium 
and recently published observations on vanadium by the alternating 
magnetic field method (Webber et al. 1949) and the ballistic method 
(Wexler and Corak 1950). 

The form of the results given above—for samples 1673 and 2088—can. 
be explained on the basis of an extension of Mendelssohn’s model of an 
“ alloy type ” superconductor (Mendelssohn 1925). This considers such 
a supe-conductor to have filaments running through it which have a much 
higher critical field (H,) than that characteristic of the bulk of the material. 
(H,). This results in the internal field of the specimen lagging behind 
that which would occur in the absence of the filaments by an amount 
(H,,) depending on the capacity of the filaments to carry resistanceless- 
currents. 

If we now also postulate a higher normal transition temperature (T,,) 
for the filaments than that for the bulk metal (T,) and an inhomogeneous 
distribution of the frozen-in flux the following results are obtained. 

The value of H, is # function of the temperature. A graph of the 
frozen-in flux, which in our use of the ballistic method is equal to 2H,, 
against temperature shows that for the two niobium specimens (nos. 1673 
and 2088) for which complete data have been obtained this function may 
be represented by H,—K(T,,—T) where K is a constant. In alternating 
magnetic field measurements the width of the transition curve is a function. 
of dH, /dT. 

The critical field of the filament material is given by the alternating 
magnetic field (and resistance) curve (curve C) and that of the bulk 
material by the higher of the ballistic curves (curve B), The transition 
temperatures and the initial slopes of the critical field curves are given. 
in Table I. 
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TABLE I, 
/ dHy dH dH 
Sample T°K. ae east om a 
ere Tort aT av dv o 
1673 8-85 9-20 1,300 1,050 45 0:55 
2088 8-35 8-75 4,000 1,900 60 0:39 


The inhomogeneous distribution of field in the material results in an: 
effective bulk demagnetizing factor (47m) for the specimen. If the higher- 
and lower of the two values of the critical field slopes given by the ballistic 
measurements (curve B and curve A) are ser ea and dH,/dT respectively 
we have dH ,/dT=(1—n) (dH,/dT). 
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Lattice Sums for Ionic Crystals. 


By F. C. FRANK, 
The H. H. Wills Physical Laboratory, University of Bristol.* 


[Received November 1, 1950. | 


THE most important of lattice sums have by now been evaluated, at 
the cost of considerable labour: but there are still many which one 
often needs or would like to know, especially for bounded lattices or for 
off-lattice points. It is not too late to weleome an improved method. 

In so far as the ionic charges can be regarded as distributed with 
spherical symmetry about the lattice points, they may be replaced (so 
far as the electrostatic field at any exterior point is concerned) by any 
other spherical distributions about the same points, or by point charges. 
This note is concerned with the evaluation of the potentials of the 
resulting lattice arrays of positive and negative point charges. 

Consider the NaCl lattice as an example. Let space be divided into 
cubical cells, each having one ion at its centre. The point charge 
(positive, say) is equivalent to the sum of two charge groups : (a) 8 point 
charges of ¢/8, placed at the corners of the cell, and (5) one charge of e 
irr 
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at the centre together with eight of —e/8 at the corners. We require the 
sum of the lattice sums for charge systems (a) and (>). The summing 
of (a) is easy, since the charges cancel at all interior points and in all 
cube faces and cube edges of the crystal. In a cuboid bounded by 
{001} faces there remain from (a) only eight charges of +e/8 at the 
corners, situated 14/3.d cube-diagonally outwards from the corner 
lattice-points, and having the signs of the corner ions. An [0, 1] step 
in the cube face has no charge, but a Kossel kink in the step has four 
charges of +e/8, situated at the corners of a square of side d (the 
interionic distance) midway between the last ion in the incomplete step 
and the first missing one ; 7. e. at the corners of the end face of the cube 
cell surrounding the last ion. The sign of these charges is that of the 
last ion. 

The charge groups (b) are equivalent to multipole systems at the 
lattice points, those of lowest order being hexadecapoles. Their potentials 
therefore decay with distance like r->. Moreover they are of alternating 
sign. The lattice sum is therefore highly convergent, and can be neglected 
save for a few lattice points close to the point of observation. However, 
it is not necessary to evaluate these multipole lattice sums as such. One 
may cut out a cuboid containing the non-negligible lattice points and 
evaluate by direct summation the potential due to the ions within it, 
together with charges of +e/8 at its corners (each opposite in sign to 
the ion in that corner), and at any other corners and kinks of the erystal. 
This is the required sum. If the calculation is to be physically realistic 
there are of course surface double layers, due to the Verwey effect, and 
additional multipoles due to ionic distortions which also have to be 
brought into account. These effects apart, the corner and kink charges 
arising from (a) reproduce exactly the charge, dipole, quadrupole and 
octopole moments of a bounded crystal. 

To test the method we use it to evaluate 


log, 2=0-693 147181 .... 


This corresponds to the potential at the next lattice point of a semi- 
infinite one-dimensional sodium chloride crystal. By our method the 
(n+-1)th approximation (requiring the addition of n+-1 terms) is 


ne 
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p=1 
n=0 (no summation) gives 0-57735, n=1 gives 0-69849, n=4 gives 
0693085, n=8 gives 0-693145 and the 6th place would then be correctly 
estimated by inspecting the alteration of terms. 

As a second test we evaluate Madelung’s sum for NaCl, which gives the 
potential at a lattice point due to the rest of the ions in an infinite crystal. 
Madelung evaluated this as 1-74(2) and Emersleben as 1-747 55(7-+2). 
Treating no ions individually we have 


2/4/38 =1:1547... 
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Treating individually the first cube shell of ions we have 
—6+12/4/2—8/1/3+2/34/3=1-748619 ... 


Taking one more shell (a total of 10 terms) gives 1-747626. The four term 
sum is sufficiently accurate for most physical applications. 

Besides possessing the merits of simple physical interpretation and easy 
application to bounded crystals or to summation for off-lattice points, the 
method appears to be more efficient for computation than any hitherto 
suggested, whether of higher or lower mathematics. The former 
methods provide admirably convergent series, but the labour comes in 
interpolating values of uncommon functions which are not tabulated 
at close intervals. The neutral shell summation methods of Evjen and 
Slater (split shells) and of Hejendahl (eccentric shells) converge much more 
slowly, having errors of nearly 2 per cent (1:77681) and 0-3 per cent 
(1-75177) after adding 9 and 12 terms respectively. 

When this method is applied to lattices of lower symmetry, the residual 
charge groups (6) will have multipoles of lower order, e. g. quadrupoles, 
making the convergence of the method slower (alternating r-* terms). 
It may then be profitable to repeat the device, so placing quadrupoles at 
kink and corner points to secure rapid convergence in the residue. 
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Some Extensions in the Use of Resistance Thermometry in the Study of 
Gaseous Explosions. 


By A. 8S. Lean, C. RountHwaitE and D. BRADLEY, 
University of Leeds*. 


[Received October 20, 1950.] 


Ix a recent communication to the Philosophical Magazine (1950, 41, 
729-30) Messrs. Daish, Fender and Woodall have raised some interesting 
points concerning a method we described in a recent paper to this Journal 


(1950, 41, 468-77) for determining flame gas temperatures by the resistance 
i emia es re oe a ee 
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thermometer method. For reasons of economy in space we were 
constrained to omit a lengthy discussion of certain aspects of the work, 
but it would now seem profitable to add a few words in amplification. 

We are entirely in agreement with them when they assert that there will 
be considerable distortion of the temperature-time curves due to lag in 
response of the wires. Indeed, it was for this reason that full allowance 
was made for the effect as they will find by referring to expression (1) on 
page 472 of our paper. In actual fact, the thermal capacities of the two 
types of wires are very little different and although the effect of the quartz 


Bigot: 
1000 
O COOLING CURVES 
Me a) PLAIN WIRE, 0:00052"D. 
Zz b) QUARTZ-COATED WIRE, 0-00052'D, 
wy 00 
a 
2 
< 400 
a 
Ww 
= ESTIMATED AIR 
2 TEMPERATURES 
= FROM 
a) PLAIN- WIRE ‘CURVE, 
b) QUARTZ-COATED WIRE CURVE. 
-200 


O 2:5 ) 75 10 12°5 IS 
TIME IN MILLISECONDS 


owing to the low conductivity of quartz, this is hardly an appreciable effect 
coating might be expected to increase the time constant of the coated wires, 
as we found from the results of a simple experiment. | 

A quartz-coated wire and then a plain platinum wire, each of 0.0005 inch 
diameter, were heated electrically up to a high temperature in stagnant air 
at atmospheric pressure within a closed vessel. A rapid-acting relay 
switch was made to cut off the heating current and switch the wire into the 
oscillograph recording circuit described in the previous paper. The 
oscillograph trace gave the cooling curves in each case and these, with the 
temperature scale, are reproduced in fig. 1. It will be noted that the 
cooling rates are not very different, in fact the time constants are 0.0100 sec. 
for the coated wire and 0.0115 sec. for the plain wire. 
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These tests, however, provide a very excel.ent check upon the accuracy 
of our method of correction, as will be seen from the broken line curves of 
fig. 1, which give the temperature of the surrounding air as deduced by 
exactly the same analysis as that employed in the estimation of flame gas 
temperatures. This temperature should have been the constant value of 
approximately 20°C. The discrepancy is not large even under the 
exacting conditions at the beginning of our records when the wire temper- 

ature was of the order of 1000°C. In fact, these experiments gave us 
full confidence in the method of correction we used. 

Thus, we disagree with the speculation of Messrs. Daish, Fender and 
Woodall that there is a large difference in the time constants of the two 
wires and that this difference could account for the catalytic heating we 
observed on the plain platinum wires, especially as full correction for the 
small difference there is in time constant was made. 

Concerning the effect of cooling of the wires by the heavier leads there is 
little we need add to what has already been said. We chose a length of 
1 inch for our wires not for ease of measurement but because repeated tests 
had shown that with lengths of 1 inch or more, the cooling effect was 
negligibly small for the quartz-coated wires. For plain wires, calculations 
and tests have shown that the mean temperature exceeds 97 per cent of 
the maximum temperature, suggesting an error of 3 per cent. But in our 
experiments this error is automatically corrected because the method of 
calibration simulates the explosive conditions very closely in so far as 
cooling by the leads is concerned (the leads, of course, do not heat up 
appreciably in 2 seconds in view of the large time constant). 

The problem of determining the true convection losses from wires in the 
explosions is one with which we have experienced difficulty. It is true 
that in the explosion there will be a certain degree of forced convection past 
the wires, as we pointed out previously. In addition, it may be that the 
equilibrium natural convection conditions have not had time to become 
fully established, although we have photographic evidence of the bodily 
convective rise of the whole flame volume in times of only a few milli- 
seconds. Nevertheless, experiments made by Mr. N. Carpenter in these 
laboratories (unpublished) have shown that the heat loss from electrically 
heated, 0.0005 inch diameter wires, in a stream of flame gas (at varying 
temperatures) hardly varies at all with gas velocity past the wires up to 
values of 30 in. persec. This phenomenon appears to indicate that the loss 
of heat is mainly dependent upon conduction through a stagnant film of 
gas around the wires. In view of this, we feel confident that our assump- 
tions regarding the convective heat loss will involve no substantial error in 
the derived temperatures and catalytic heating rates. 
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CXVI. Notices of New Books and Periodicals received. 


Structure of Molecules and the Chemical Bond. By J. K. Syrxry and 
M. E. Dyarxkra, translated and revised by M. A. Partridge and D. O. Jordan. 
(Butterworths Scientific Publications.) 


T's book describes the properties of the chemical bond and interprets them in 
terms of quantum mechanics. 

It is clearly addressed mainly to experimental workers and therefore starts 
with an account of the quantum mechanics necessary to interpret the experi- 
mental results, the less essential part of the mathematics being relegated to an 
appendix. There are two different approximate solutions of the wave equation 
that will give a semi-quantitative interpretation of the chemical bond; these 
are the electron-pair bond plus resonance approximation, and the molecular 
orbital approximation. Both these approximations are described and compared, 
but in the subsequent discussion of experimental results it is usually the first 
approximation that is used. The latter chapters of the book consist of a large 
amount of experimental data on the chemical bond and the interpretation of 
this data in terms of the theory given in the first chapters. The subjects con- 
sidered include bond lengths, bond energies, dipole moments, van der Waals 
forces, and the general arrangement of bonds in complex molecules and crystals. 
A minor criticism is that in Chapter 10 the word polarization and the symbol P 
are used in an unusual sense that makes the discussion rather puzzling at first. 

Although this book is addressed mainly to experimentalists, the theoretical 
worker will also find the collection of experimental data of great value. A. F. D. 


Calculating Instruments and Machines. By D. R. Hartren. (Urbana> 
University of Illinois Press.) 


THE English edition of this book is now obtainable, price 21s. Published by 
Cambridge University Press, London. 


BOOKS RECEIVED. 


A new scientific publication has been received— Geochimica et Cosmochimica 
Acta’, edited with the assistance of an Advisory Board. Published by 
Butterworth-Springer Ltd. (London). Price per vol. (approx. 300 pp.) 
£3 10s. Od. 


[The Editors do not hold themselves responsible for the views. 
expressed by their correspondents. | 
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Apams (N.), temporary increase 
in neutron compound of cosmic 
rays, 503. 

—, (with Braddick, H. J. J.), 
505. 

Atkins (K. R.) and Osborne (D. V.), 


the velocity of second sound 
below 1° K., 1078. 
Alpha particles, short range 


emission.—II. and III. (Burcham 
and Freeman), 337, 873. 

; IV. (Freeman), 1225. 

—,, , from break up of &Li 

(Baxter, Burcham and Paul), 937. 

Alternant hydrocarbons, molecular- 


orbital calculations (de Heer), 
370. 

Atmospheric absorption of heat 
radiation by water vapour 


(Cowling), 109. 

Atmospherics, use to study very 
long radio waves (Gardner), 1259. 

Atomic oxygen in CO explosions 
(Leah, Rounthwaite and Bradley), 
478. 

Avery (D. G.), some optical pro- 
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